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Abstract 


■  i 

Application  of  optical  signal  processing  techniques  to  adaptive 
array  processing  was  investigated.  The  common  adaptive  array  algorithms 
were  reviewed.  Adaptive  array  processing  generally  requires  some  form 
of  correlation  and  the  characteristics  of  optical  correlators  were 
compared  leading  to  the  selection  of  the  time  integrating  optical 
correlator  for  further  investigation.  An  adaptive  array  processor 
using  time  integrating  optical  correlators^  which  was  first  proposed  by 
David  Casasent,  was  analyzed  to  determine  applicability  of  this  signal 
processing  technique  to  a  general  adaptive  array  environment.  This 
processor  was  found  to  be  similar  to  a  comparable  maximum  likelihood 
approach.  Limitations  on  the  performance  of  this  processor  in  terms 
of  ambiguity  and  resolution  were  found  and  these  results  were  con> 
firmed  by  computer  simulation  of  the  adaptive  array  environment,  the 
processor,  and  the  resulting  antenna  response. 


I.  Introduction 


An  adaptive  array  rejects  interference  by  adjusting 
the  amplitude  and  phase  response  at  each  of  its  elements  to 
place  nulls  In  the  antenna  pattern  in  the  directions  of  the 
interference.  Adaptive  processing  becomes  complex  as  the  number 
of  adaptive  weights  increases.  The  use  of  optical  signal  processing 
offers  numerous  advantages  over  digital  or  electronic  analog  pro¬ 
cessors  such  as  speed  and  parallel  processing. 

This  thesis  presents  the  results  of  an  investigation  of 
the  application  of  optical  signal  processing  to  adaptive  arrays. 

This  investigation  began  with  very  few  constraints  on  the  directions 
to  be  pursued  so  a  considerable  effort  was  necessary  to  determine 
which  areas  of  the  problem  deserved  the  most  attention.  The  nature 
of  the  adaptive  array  probiem  was  reviewed  to  provide  an  understanding 
of  its  signal  processing  needs.  Optical  signal  processing  was  then 
surveyed  to  determine  the  possibilities  of  meeting  these  processing 
needs.  Chapter  II  of  this  thesis  reviews  the  adaptive  array  back¬ 
ground  and  discusses  some  of  the  algorithms  which  have  been  developed 
for  adaptive  processing.  A  common  factor  in  many  of  the  adaptive 
algorithms  is  correlation  and  Chapter  III  provides  an  overview  of 
the  diverse  possibilities  for  optical  correlators.  The  time  inte¬ 
grating  form  of  optical  correlator  was  chosen  as  the  best  for  the 
adaptive  array  problem  and  Chapter  IV  presents  the  characteristics  of 


each  part  of  this  correlator  with  «r.f«hasis  on  the  performance  iimiting 
aspects  of  each.  An  adaptive  array  processor  utilizing  time  integra¬ 
ting  optical  correlation  which  was  first  presented  by  Casasent  (Ref  1,2) 
and  variations  on  this  processor  are  the  subjects  of  Chapter  V.  The 
concepts  upon  which  this  processor  is  based  is  to  use  a  correlation 
of  the  signal  present  at  the  elements  of  the  adaptive  array  to 
estimate  the  directions  of  the  interference  sources.  This  process 
is  comparable  to  a  maximum  likelihood  estimate  of  these  directions. 

These  directions  must  be  used  to  determine  element  weights  to  null 
these  direction.  An  efficient  approach  to  this  post  processing  task 
is  presented  which  has  considerable  flexibility  in  terms  of  imple¬ 
mentation.  As  part  of  this  thesis  the  time  integrating  optical 
correlator  and  associated  post  processing  were  simulated.  Details 
of  this  simulation  are  included  in  Chapter  VJ  and  the  Appendix,  and  the 
results  of  this  simulation  and  conclusions  concerning  optical  signal 
processing  for  adaptive  arrays  and  the  performance  of  this  particular 
time  integrating  processor  are  contained  in  Chapter  VI. 


II.  The  Adaptive  Array 


Introduction 

An  adaptive  array  (AA)  rejects  interference  by  placing 
nulls  in  the  receiving  antenna  pattern  in  the  directions  of  inter¬ 
fering  signals.  Ideally,  this  is  accomplished  without  losing  sen¬ 
sitivity  in  the  direction  of  the  desired  signal.  The  adaptive 
array  compensates  automatically  for  spatial,  temporal,  and  spectral 
variations  in  the  signal  plus  noise  environment.  (Ref  3)2).  Optimum 
matching  of  the  temporal  response  of  the  system  to  the  signal  and 
noise  waveforms  contributes  to  the  depth  of  the  spatial  pattern  nulls 
(Ref  3*3 )•  The  temporal  aspect  of  the  system  involves  processing 
the  incoming  signals  and  continually  adjusting  the  pattern  tqthe 
changing  environment.  The  active  part  of  the  adaptive  array  is  the 
processor  which  uses  some  a  priori  information  about  the  signal 
or  noise  and  the  signal  plus  noise  inputs  from  the  array  to  deter¬ 
mine  the  amplitude  and  phase  weighting  of  each  element  of  the  array 
which  is  optimum  according  to  some  given  criterion.  The  criterion 
chosen  and  the  constraints  which  are  placed  on  the  optimization  lead 
to  different  forms  of  the  adaptive  array  processor. 

This  chapter  will  first  present  some  background  including 
the  notation,  assumptions,  and  basic  characteristics  upon  which 
later  discussions  of  the  adaptive  array  will  be  based.  The  re¬ 
mainder  of  the  chapter  will  present  an  overview  of  the  common 
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algorithms  used  for  AA  processing  and  the  variations  for  different 
a  priori  knowledge  as  well  as  discussing  important  performance 


measures. 

General  Structure 

The  general  form  of  the  adaptive  array  as  shown  in  Fig.  1 
consists  of  the  antenna  array,  fixed  processing,  element  weights, 
and  the  weight  controlling  system  or  AA  processor.  A  given  system 
may  combine  these  in  different  ways  and  different  processors  may  use 
some  subset  of  the  inputs  shown.  The  notation  which  will  be  used 
in  this  thesis  for  AA  discussions  is  given  in  Table  1. 

Array  Configuration 

The  antenna  elements  may  be  arranged  lineaeVy,;  circularly, 
or  even  randomly.  The  performance  of  the  array  can  be  limited  or 
enhanced  by  different  geometries  (Ref  3:23)  and  the  form  of  the 
processor  is  different  for  different  geometries.  A  linear  array  of 
evenly  spaced  elements  was  chosen  for  this  thesis.  This  choice 
leads  to  a  simpler  processor  for  the  processing  approach  investigated 
Variations  in  element  spacing  for  the  linear  array  result  in  varying 
null  widths  and  thus  are  important  to  AA  performance.  (Ref  3:20). 
However,  the  form  of  the  processor  is  not  affected  so  no  further 
assumptions  will  be  made  on  this  point. 

Antenna  Response 

The  linear  array  with  the  equally  spaced  elements  exhibits 
an  equal  dealy  from  element  to  element  for  a  signal  for  a  directional 


Adaptive  Array  Processing  Structure 


source  which  is  sufficiently  distant  to  be  represented  as  plane  wave 
fronts.  The  delay  for  a  source  at  an  angle  8  off  of  boresight  is 
given  by 

T  -  -  sine  (1) 

c 

where  d  is  the  element  to  element  spacing  and  c  is  the  speed  of  sig-> 
nai  propagation.  This  equation  is  easily  seen  from  Fig.  2. 


Fig  2.  Eiement-to-Eiement  Time  Delay 
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Table  1.  Adaptive  Array  Notation 

Number  of  elements  1 

A  given  array  element  (element  zero  is  far  right) 

Number  of  directional  noise  sources 

A  given  noise  source 

Center  to  center  separation  of  elements 

Wavelength  of  radiation 

Signal  received  at  element  n 

Angle  of  source  k  from  bore  sight. 

Angle  of  source  off  of  bore  sight 
Delay  of  source  k  from  element  to  element 
A  vector  of  the  complex  weights  applied  to  each  element 
The  weight  applied  to  the  nth  element 
Signal  process  from  noise  source  k  as  received  at  ele¬ 
ment  zero 

Desired  signal  as  received  at  element  zero 
Signal-plus  noise  covariance  matrix 


For  a  source  with  frequency,  f,  the  phase  delay  is  thus 


(j>  -  2it  f  i  sine  (2) 

c 

The  AA  multiplies  the  signal  at  element  n  by  some  weight,  so  the 

far  field  response  of  the  array  can  be  found  to  be 

N 

B(0)  m  I  u  exp(jn2irf  sln0)  (3) 

n-0  " 

where  far  field  can  be  taken  mean  distances  so  much  greater  than  the 
largest  dimension  of  the  array. 

Input  and  the  Covariance  Matrix 

The  general  form  of  the  input  and  the  covariance  matrix 
are  presented  here  for  later  reference.  The  input  signal  r^(t) 
element  n  can  be  written  in  the  following  form: 

K 

rn(t)  -  Sj(t-nTj)  +  S  S,^(t-nT|^)  (k) 

k»1 

for  a  desired  signal  at  some  angle  0^  from  bore  sight  and  K  directional 
noise  sources  at  angles  0|^.  Note  that  S^(t)  and  S|^(t)  are  the  sig* 
nals  as  received  at  element  zero,  the  far  right  element  looking 
towards  bore  sight  and  that  they  may  be  deterministic  or  random  in 
nature.  The  delay  of  each  signal  at  element  n  is  thus  nr  where  t  is 
as  defined  in  eq.  (it). 

The  covariance  matrix,  N,  is  made  up  of  the  covariance 


terms 


The  resulting  matrix  is  fundamental  to  AA  theory.  The  terms  are 
determined  by  the  power  and  spatial  distribution  of  the  noise  sources 
The  desired  weight  vector,  U,  can  be  shown  to  be  the  result  of  a 
vector  matrix  product, 

W  -  ^  (6) 

where  ^  is  a  steering  vector  which  determines  the  direction  of  the 
main  beam  (Ref  4:586  and  1:6). 

Algorithms 

This  section  will  brief ly  review  the  common  algorithms  which 
have  been  developed  for  AA  processing.  Most  of  the  algorithms  which 
have  been  developed  to  find  the  weight  vector,  W,  for  a  given  noise 
environment  can  be  derived  as  an  optimal  solution  based  on  a  parti¬ 
cular  performance  criterion.  These  algorithms  and  their  associated 
performance  criteria  are: 

LHS  -  Minimizes  the  mean  square  error  between  the  array 
output  signal  and  a  reference  signal 

MSN  -  Maximizes  the  output  signal  to  noise  ratio 

ML-  Finds  a  maximum  liklihood  estimate  of  the  signal 
information 


MNV  -  Minimizes  the  output  noise  variance 


Other  algorithms  which  have  been  developed  include  the  SMI  or  sample 

e 

matrix  inversion  approach  and  the  RS  or  random  search  method. 

IMS,  The  LMS  algorithm  was  developed  by  B.  Widrow  (Ref  5)  and  it 
has  probably  received  the  most  widespread  application  (Ref  3:9)« 

There  are  several  variations,  but  in  general,  a  reference  signal  is 
subtracted  from  the  output  and  the  result  is  correlated  with  the 
incoming  signal  at  each  element.  Many  references  are  available 
detailing  the  theory  and  performance  of  this  algorithm  (Ref  5,6,7). 

MSN.  The  maximum  signal  to  noise  ratio  (MSN)  algorithm  was  developed 
by  P.U.  Newells  and  S.P.  Applebaum  (Ref  8).  This  algorithm  is  the 
result  of  radar  side  lobe  cancelling  work.  A  generalized  output 
signal  to  noise  ratio  is  maximized  by  means  of  a  steering  signal 
at  the  element  weights.  Both  the  MSN  and  LMS  algorithms  derive 
tHeir  weight  control  from  a  correlation  between  the  signals  at  each 
of  the  elements  and  both  converge  toward  the  same  solution  (Ref  4). 

ML.  The  maximum  likelihood  (ML)  algorithm  produces  an  estimate 
of  the  desired  signal  which  is  the  most  likely  waveform.  This 
optimization  requires  known  second  order  statistics  of  the  noise 
and  generally  assumes  a  Gaussian  distribution  (Ref  6:3-4,  3-10). 

This  algorithm  requires  no  knowledge  of  the  signal  except  that  it  is 
deterministic. 

MNV.  The  minimum  noise  variance  (MNV)  algorithm  minimizes  the  out¬ 
put  noise  variance  (Ref  6:3-2).  This  algorithm  assumes  known 
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direction  of  arrivai  and  zero  mean  noise  sources  such  that  the  out¬ 
put  after  correction  for  the  direction  of  arrival  is  an  unbiased 
estimate  of  the  desired  signal.  This  method  also  requires  no  know¬ 
ledge  of  the  signal  and  only  second  order  statistics  of  the  noise. 

SMI .  The  sample  matrix  inversion  (SHI)  algorithm  is  an  example  of 
a  higher  order  processing  algorithm.  A  direct  calculation  of  the 
optimum  weights  is  made  based  on  sanpling  the  signal  at  each  element 
and  inverting  the  resulting  noise  covariance  matrix  (Ref  3:11)*  Thi 
method  produces  the  weights  directly  but  Is  computations I y  complex. 

RS.  The  random  search  (RS)  algorithm  varies  the  weights  in  a  random 
or  nearly  random  manner  in  such  a  way  as  to  continually  improve  some 
performance  measure  (Ref  3:11  and  Ref  6:4-2).  This  is  basically  a 
refined  trial  and  error  approach.  This  has  the  advantage  of  sim¬ 
plicity,  however,  some  guidance  must  be  added  to  this  search  and 
in  general,  the  search  techniques  do  not  lead  to  practical  circuitry 
(Ref  6:4-5). 

Optimization  Techniques  (Ref  6),  The  technique  employed  by  each 
algorithm  to  optimize  the  array  weights  directly  affects  the  rate  of 
convergence.  These  algorithms  can  basically  be  categorized  as  search 
(zero  order),  gradient  (first  order),  or  higher  order  techniques. 

Most  algorithms  which  have  been  investigated  Including  the  LMS, 

MSN,  ML,  and  MNV  use  the  first  derivative  of  the  performance  cri¬ 
terion  and  converge  on  the  optimum  point  by  adjusting  the  weights 
to  move  In  the  direction  of  steepest  descent  towards  the  optimum 


weights.  The  rate  of  convergence  is  determined  by  the  size  of  the  ad¬ 
justment  made  at  each  step.  The  size  of  this  adjustment  is  the  gain 
of  the  feedback  ioop  and  as  with  any  controi  loop  the  choice  of  this 
parameter  is  critical  to  the  speed  or  adaptation  time  and  the  stabiiity 
of  the  processor.  Higher  order  methods  effectively  use  some  informa¬ 
tion  about  the  second  or  higher  derivatives  and  thus»  these  methods 
are  able  to  produce  a  set  of  weights  which  are  optimum  by  performing 
only  one  iteration.  The  trade-off  which  must  be  considered  is 
between  complexity  of  a  higher  order  systems  and  slow  convergence  of 
iower  order  systems. 


Summary  of  Algorithms,  The  aigorithms  discussed  provide  a  generai 
overview  of  the  types  of  AA  processors.  There  are  untoid  variations 
on  these  approaches  and  others  which  were  not  discussed.  Ai though 
the  criteria  which  these  algorithms  optimize  are  different,  they 
ail  approach  the  same  resuits.  One  of  the  major  differences  is 
difficulty  of  implementation  and  another  is  the  rate  of  convergence 
to  the  optimum  set  of  weights.  Convergence  has  been  discussed  here 
and  the  reader  is  referred  to  Ref,  6,  Chapter  X  for  a  discussion  of 
implementation. 


A  Priori  Knowledge 

The  form  of  the  AA  processor  is  dependent  upon  the  informa¬ 
tion  which  is  known  a  priori.  The  AA  basically  eliminates  unwanted 
signals  by  adding  two  or  more  signals  which  have  phase  differences 
such  that  their  sum  is  zero.  To  discriminate  between  signal  and 
noise  some  a  priori  knowledge  must  be  available  (Ref.  3:13).  This 
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Information  might  be  the  waveform*  spectral  information*  directions 
of  arrival  (DOA)  *  polarizations*  or  power  levels.  The  knowledge 
might  be  about  the  noise  or  the  signal  and  of  course*  many  combinations 
are  possible. 

Performance  Measures  (Ref,  3i17“19) 

Adaptive  arrays  differ  greatly  from  one  algorithm*  or 
implementation  or  set  of  a  priori  information  to  another.  As  a  re¬ 
sult*  the  performance  measures  used  to  characterize  them  are  diverse 
and  comparison  can  be  difficult.  Each  adaptive  array  emphasizes  cer¬ 
tain  criteria  which  are  important  to  a  particular  application.  Some 
of  these  characteristics  include  the  bandwidth  over  which  the  nulls 
exist*  the  number  of  array  elements*  signal  optimization*  or  all  sig¬ 
nal  suppression*  adaptation  time*  and  number  of  sources  to  be  nulled. 
Operating  characteristics  such  as  the  ability  to  retain  nulls  for 
intermittant  Jammers*  called  null  memory  and  the  methods  used  to 
acquire  the  desired  signal  are  also  important. 

A  common  point  for  comparison  is  null  depth.  This  measure 
is  analogous  to  processing  gain  in  a  spread  spectrum  system.  Direct 
comparison*  however*  of  this  measure  may  be  inappropriate  due  to 
differences  in  bandwidth  or  adaptation  time.  A  useful  measure  of  an 
adaptive  array  is  the  improvement  of  signal  to  noise  ratio  at  the 
output  since  more  factors  are  included  in  its  determination,  in 
general*  the  specific  needs  of  the  application  and  a  number  of  the 
performance  measures  which  have  been  mentioned  must  be  taken  into 
consideration  for  any  valid  comparison. 


Summary 

This  chapter  has  presented  some  background  information  con¬ 
cerning  adaptive  arrays.  The  general  form  of  the  array  and  its  pro¬ 
cessor  were  presented  and  some  related  notation  and  the  response 
at  the  array  output  and  input  were  introduced.  A  number  of  pro¬ 
cessing  algorithms  were  briefly  discussed  and  some  related  topics 
were  included.  Finally,  an  overview  was  made  of  some  of  the  per¬ 
formance  measures  used  for  adaptive  arrays.  This  background  will 
be  the  basis  upon  which  presentation  of  the  time  integrating  optical 
correlator  for  adaptive  processing  will  be  made. 
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IN.  Optical  Signal  Processing 


Introduction 

Optical  signal  processing  (OSP)  Is  one  of  the  answers  to  the 
technological  problem  of  processing  large  quantities  of  signal  informa¬ 
tion  for  real  time  applications.  Optical  signal  processing  architectures 
and  algorithms  begin  with  electrical  signals,  convert  the  electrical 
signals  to  light  signals,  process  the  light,  and  convert  it  back  to 
the  electrical  signal  which  is  the  desired  product.  This  chapter  will 
present  some  of  the  background  needed  for  a  good  understanding  of  OSP 
and  then  some  of  the  architectures  which  have  been  developed  will  be 
briefly  presented.  The  emphasis  here  will  be  entirely  on  correlating 
processors  since  correlation  is  a  common  factor  in  most  of  the  adaptive 
array  algorithms. 

Background 

Light  Propagation.  This  section  will  briefly  discuss  some  of  the  light 
properties  which  are  used  in  OSP.  The  propagation  of  light  as  a  wave 
can  be  represented  by  the  Huygen-Fresnel  equation.  This  equation  can 
be  simplified  in  the  far  field  or  in  the  focal  plane  of  a  lens  and 
the  result  is  the  Fraunhoffer  equation.  In  the  focal  plane  of  a  lens 
the  Fraunhoffer  equation  reduces  to  a  constant  times  the  two  dimen¬ 
sional  Fourier  transform  of  the  amplitude  distribution  of  the  light 
entering  the  lens.  This  is,  perhaps,  the  most  important  optical 
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property  used  in  OSP.  Using  this  property  a  transformation  can  be 
made  from  the  time  domain  to  the  frequency  domain  at  the  speed  of 
light  with  only  a  simple  lens  system.  This  same  property,  however, 
aiso  limits  the  ultimate  resoiution  obtainabie  with  any  common  opti¬ 
cal  system.  The  smallest  resoluable  area  is  limited  by  diffraction 
through  the  optical  system.  The  resolution  iimit  for  an  aperture 
whose  size  is  given  by  d  in  a  particular  direction  is  about 

6-^  (7) 

in  that  same  direction  where  is  the  wavelength  of  light  being  used 
and  L  is  the  distance  from  the  aperture.  The  properties  discussed 
to  this  point  can  all  be  obtained  through  Fourier  optics  (Ref  9). 

Modulators 

All  OSP  systems  can  be  broken  down  Into  three  parts; 
the  light  source,  a  light  modulator  or  modulators,  and  the  detectors 
and  the  optics.  These  will  all  be  discussed  more  later  but  some  of 
the  common  light  modulators  will  be  presented  here  in  order  that 
later  discussions  may  be  limited  to  only  one  type.  Some  of  the  de¬ 
vices  which  have  been  used  are  listed  here: 

1.  Acousto-optic  cells 

2.  Liquid  crystal  light  value  (Ref.  10). 

3.  Thermo  plastic  optical  phase  devices  (Ref.  10) 

4.  Electro-optic  cells 

Only  the  first  three  have  the  ability  to  spatially  modulate  light 
and  only  the  acousto-optic  (AO)  cell  has  the  speed  which  is  needed 
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for  signal  processing  applications.  Although  technoiogical  improve* 
ment  may  make  other  systems  viable  in  the  future,  oniy  acousto-optic 
systems  will  be  investigated  in  this  thesis. 


Correlator  Characteristics.  The  following  sections  will  present  some 
of  the  OSP  architectures  which  have  been  proposed  to  perform  the 
operation  of  correlating  two  signals.  The  correlation  to  be  performed 
is  given  here 

S(t)  -J  S,(f)  S2(t+t')  df  (8) 

o 

Commonly,  one  of  these  signals  Is  a  reference  signal  and  the  other  is 

a  received  signal  containing  both  signal  and  noise.  Some  basic 

characteristics  of  this  correiation  are  important  in  understanding  a 

comparison  of  the  processors  to  be  discussed.  The  first  is  the  width 

of  the  correiation  peak  output. 

The  correlation  of  a  simple  pulse  of  length  T^,  with  a 

similar  pulse  produces  a  correlation  of  width  2T.  This  signal  has  a 
2 

sine  frequency  distribution  with  a  main  lobe  width  of  2/T^  as  shown 
in  Fig.  3*  The  width  of  the  correlation  peak  is  the  reciprocal  of  the 
main  lobe  bandwidth  yielding  an  approximation  for  the  relationship 
between  bandwidth  and  correlation  width  of 

W  4  i  (9) 

where  B  is  the  bandwidth  of  the  signals  being  correlated  (Ref  11:138). 

Another  important  characteristic  is  the  power  signal-to-noise 
ratio  of  the  correlation  output  as  compared  to  the  signal-to-noise  of 
the  noisy  signal  S^(t) 
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SNR  *  BT  SNR 


1 


(10) 


where  SNRj  is  the  signal-to-noise  of  S^(t)»  T  is  the  correiation  time, 
and  BT  is  the  time  bandwidth  product  of  the  correlator  (Ref.  11:137)* 

The  gain  is  reduced  if  the  bandwidth  of  the  signal  is  less  than  the 
bandwidth  of  the  correlator  or  if  the  duration  of  the  signal  is  shorter 
than  the  integration  time.  The  correlator  thus  exhibits  a  gain  in  signal- 
to-nolse  called  correlation  gain  or  processing  gain,  equal  to  the  time 
bandwidth  product  of  the  processor.  The  correlation  time,  T,  is  called 
the  time  window.  Another  parameter  called  the  range  window  is  the  time 
error  allowed  between  the  received  and  reference  signal.  The  received 
signal  may  in  fact  be  composed  of  the  reference  signal  delayed  by  some 
quantity,  T,  plus  the  noise  signal 

S^(t)  -  S2(t-T)  +  n(t)  (11) 
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If  this  error,  t,  is  greater  than  the  time  window,  T,  no  overlap  of  the 
signals  will  occur  where  the  reference  portion  of  the  received  signal 
and  the  reference  signal  coincide. 

These  four  parameters 

1.  Bandwidth 

2.  Time  with  product 

3.  Time  window 

4.  Range  window 

can  be  used  to  compare  the  OSP  configurations  which  are  available  for 
AA  processing. 

OSP  Architectures 

Spatial  Integrating.  The  simplest  and  the  tranditional  form  of  acousto¬ 
optic  correlator  performs  the  integration  operation  over  a  spatial 
coordinate.  This  correlator  first  produces  a  light  wave  whose  ampli¬ 
tude,  phase  or  intensity  is  a  spatial  representation  of  one  of  the 
signals  to  be  correlated.  This  light  wave  is  then  further  modulated 
by  the  second  signal  which  produces  the  product  along  a  spatial  coordi¬ 
nate.  Integration  along  this  coordinate  is  performed  by  focusing  this 
product  onto  a  point  detector.  The  spatial  light  signals  can  be 
produced  by  AO  modulators  and  as  the  signals  propagate  through  the 
modulators,  the  signal  from  the  detector  is  the  correlation  for 

different  time  delays.  A  schematic  of  this  type  of  correlator  is 

4 

shown  in  Fig.  4,  Notice  that  the  signal  into  the  second  cell  must 
be  the  reference  signal  time  reversed  in  order  for  the  correlation 
to  be  formed.  The  two  AO  cells  can  in  fact  be  combined  into  one 
cell  to  produce  a  very  compact  unit  which  requires  very  few  optical 
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elements  and  is  easily. aligned  (Ref.  11:140).  Many  variations  on 
this  general  idea  are  in  fact  possible.  For  a  known  reference  signal 
a  fixed  mask  may  replace  one  of  the  AO  ceils.  It  should  be  pointed 
out  that  the  figure  shown  is  only  a  schematic.  Variations  which  may 
be  made,  involve  such  things  as  using  coherent  or  incoherent  light. 
With  coherent  light,  a  complex  correlation  may  be  performed  using 
heterodyne  detection  (Ref.  12:70).  For  this  type  of  detection,  the 
modulators  actually  phase  modulate  the  light  wave  and  this  phase 
modulation  is  converted  to  amplitude  modulation  (Ref.  12:67). 

Other  spatial  integrating  architectures  use  the  Fourier  transform 
property  of  light  diffraction.  A  unique  form  of  the  spatial  corre¬ 
lator  called  a  chirp  correlator  detects  the  time  of  arrival  of  linear 


20 


frequency  modulated  pulses  (chirps).  The  acoustic  wave  modulated 
by  the  chirp  acts  as  a  Fresnel  zone  plate  which  focuses  light  on  a 
detector  providing  the  desired  correlation  (Ref.  11:14). 

The  system  specifications  for  these  different  forms  of 
spatial  correlators  are  largely  determined  by  the  AO  cell.  The  maximum 
signal  bandwidth  is  the  bandwidth  over  which  the  AO  ceil  operates, 
efficiently  and  with  acceptable  distortion.  Detectors  are  available 
which  do  not  limit  this  bandwidth.  The  time  window  is  limited  to  the 
time  window  of  the  AO  cell  since  integration  is  performed  spatially  along 
the  length  of  the  cell.  The  time  bandwidth  product  or  gain  of  the 
correlator  is  the  product  of  these  and  thus  is  determined  solely  by 
the  AO  ceil.  AO  ceils  can  be  produced  with  time  bandwidth  specifications 
as  high  as  1600  (Ref.  13:2).  The  range  window  of  the  correlator  as 
shown  in  Fig.  4  is  limited  to  the  time  window  of  the  AO  cell.  This 
is  because  even  for  very  short  pulses  the  signals  must  be  present  in 
both  cells  to  produce  a  correlation.  The  second  cell  may  be  replaced 
by  a  reference  mask  for  a  known  signal  and  the  resulting  range  window 
is  infinite.  The  range  window  for  the  chirp  correlator  is  also  in¬ 
finite  since  no  reference  is  used. 

Time  Integrating.  The  basic  form  of  the  time  integrating  acousto¬ 
optic  correlator  is  shown  in  Fig.  5.  The  intensity  of  the  light 
source  is  modulated  by  S^(t)  and  the  product.  S|(t)  S2(t'«-t')  is 
formed  In  .  the  AO  ceil.  The  time  delay  t'  is  determined  by  the 
position  in  the  ceil  and  velocity  of  propagation  of  the  acoustic  wave. 

The  product  is  then  integrated  in  time  at  the  detector  array  which  is 
in  general  a  CCD  array. 
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Fig.  5.  Time  integrating  Correlator 


To  better  understand  the  processes  taking  place  In  the 
time  integrating  processor,  it  is  helpful  to  trace  the  signal  pro¬ 
gression  through  the  system  as  shown  in  Fig.  6.  The  system  being 
described  uses  a  modulated  point  source  collimated  by  a  cylindrical 
lens.  Shown  in  part  'a'  of  that  figure  are  the  light  pulses  traveling 
towards ithe  AO  cell  and  another  pulse  which  is  the  input  to  the  AO  cell 
In  'b*  the  light  pulses  have  Just  entered  the  cell.  In  'b'  and  'c'. 
the  light  pulse  which  is  the  result  of  the  second  signal  is  shown 
exiting  the  cell  along  with  the  first  pulse.  The  second  signal  is 
drawn  slightly  taller  for  identification  and  the  product  formed  is 
shown  in  part  'c*  by  the  cross  hatched  area.  The  detectors  integrate 
the  product  over  time  and  output  the  result  at  the  end  of  the  integra¬ 
tion  time.  Shown  in  part  'd*  at  the  detectors  is  the  output  for  the 
case  shown.  Due  to  the  delay  as  the  pulse  travels  through  the  AO  cell 
the  pulses  are  synchronized  in  time  only  along  one  ray  of  light,  thus 
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Time  Integrating  Correlator  Operation 


producing  the  peak  output  at  this  point. 

Another  time  integrating  optical  correlator  (TIOC)  sche¬ 
matic  is  shown  in  Fig.  7.  This  architecture  uses  two  AO  ceils  (or 
one  cell  with  two  acoustic  channels).  The  product  is  formed  as  the 
light  is  modulated  by  each  of  the  cells  and  the  result  is  integrated 
in  time  by  the  detector  array. 

The  maximum  signal  bandwidth  for  the  TIOC  is  also  limited 
to  the  bandwidth  of  the  AO  cell.  The  time  window  of  this  processor  is 
the  time  over  which  the  detector  is  allowed  to  integrate.  This  time 
window  is  thus  infinite,  however,  the  time  bandwidth  product  or  gain 
in  signal  to  noise  ratio  is  limited  by  detector  noise.  The  maximum 
signal  to  noise  voltage  ratio  for  the  CCD  array  is  given  by  its  dynamic 
range,  D.  The  time  bandwidth  product  is  effectively  D  (Ref,  11:143), 

In  fact,  there  are. other  noise  and  bias  sources  which  reduce  the  gain 
but  detectors  are  available  with  dynamic  ranges  of  10^  (10^  with  cooling) 
so  very  high  gains  are  possible.  The  range  window  of  the  TIOC  is 
limited  to  the  time  window  of  the  AO  cell. 


Fig.  7.  Two  Cell  Time  Integrating  Correlator 


To  determine  the  ability  of  these  correlators  to  provide 
AA  processing  needs,  a  comparison  of  their  basic  correlation  para¬ 
meters  can  be  made.  A  summary  of  correlator  specifications  is  given 
here  in  Table  2  as  taken  from  (Ref.  11). 

Table  2.  Correlator  Comparison 

Time 


Spatial  integrating 

Band¬ 

width 

Time 

Window 

Bandwidth 

Product 

Range 

Window 

a.  Fixed  reference  mask 

B 

T 

Bt 

GO 

b.  Two  cell  correlator 

B 

T* 

Bt*- 

T 

c.  Chirp  correlator 

B 

T 

Bt 

00 

Time  integrating 

B 

00 

T 

B  a  bandwidth  of  AO  cell 
”■  X  "  tJ'"®  window  of  AO  cell 
D  ■  detector  dynamic  range 


^Assuming  complete  overlap  of  the  two  signals  for  short  pulses 

The  space  integrating  and  time  integrating  systems  each  have  advantages 
and  disadvantages.  The  space  integrating  system  has  a  limited  time 
bandwidth  capability  but  is  capable  of  searching  over  a  large  range 
window  or  delay  in  time  of  arrival.  The  time  integrating  system  allows 
processing  of  large  time  bandwidth  product  signals  thus  producing 
significant  gain;  however,  the  range  window  is  limited  (Ref.  14:42). 

For  the  AA  processing  needs,  the  signal  delay  Is  never  longer  than  the 
propagation  delay  between  the  most  widely  separated  elements.  Time 


windows  for  AO  cells  are  commonly  several  microseconds  so  the  range 
window  is  not  an  important  consideration  for  AA  processing.  The  time 
integrating  optical  correlator  shows  the  greatest  potential  for  AA 
processing  particularly  for  signal  environment,  with  large  time  band¬ 
width  signals.  The  time  integrating  systems  achieve  processing  gain 
over  the  full  duration  of  a  signal  (up  to  the  limits  set  by  dynamic 
range)  whereas  the  space  integrating  system  is  limited  to  durations  of 
several  microseconds.  Because  of  these  observations,  further  investigation 
of  OSP  in  this  thesis  is  limited  to  acousto-optic  time  integrating  cor¬ 


relators 


IV.  Time  integrating  Optical  Correlators 

The  basic  form  and  the  correlation  specifications  for  the  time 
integrating  optical  correlator  (TiOC)  were  presented  in  the  last  chapter. 
This  chapter  investigates  the  characteristics  of  each  part  of  the  TIOC 
with  particular  emphasis  to  those  aspects  such  as  noise  which  limit 
the  overall  performance. 

The  parts  of  the  TiOC  which  will  be  examined  are  the  source,  the 
AO  modulator,  and  the  detector.  Other  parts  of  the  TIOC  which  will  not 
be  discussed  are  the  optics  and  the  electronics  which  drive  the  AO  cell 
and  process  the  detected  signals.  The  lenses  required  for  these  sys¬ 
tems  may  be  obtained  with  sufficient  quality  that  they  do  not,  in 
general,  limit  the  performance  of  the  system.  It  should  be  noted, 
however,  that  very  high  quality  optics  must  be  specified  for  heterodyne 
detection  where  phase  information  must  be  maintained. 

AO  Modulators 

Theory.  This  section  will  present  some  of  the  basic  theory  of  acousto¬ 
optic  cells  necessary  to  understand  the  operating  limitations  and  first 
order  design  criteria  of  these  cells.  Bragg  diffraction  can  be  thought 
of  in  somewhat  the  same  way  as  X-ray  diffraction  in  a  crystal  with  the 
atomic  planes  replaced  by  the  planes  of  compression  and  rarefaction 
introduced  by  ultrasonic  sound  waves  (Ref.  15:4).  The  velocity  of 
the  sound  waves  causes  the  diffracted  light  to  be  Doppler  shifted.  As 


shown  in  Fig.  8,  a  beam  of  light  is  directed  into  the  path  of  the 
acoustic  wave  inside  the  ceil.  At  certain  angles  of  incidence  +  oig,  a 


,  Fig,  8.  Bragg  Cell  Diffraction 

diffracted  beam  is  produced  whose  direction  differs  by  The  angle 

Og  is  called  the  Bragg  angle  and  is  given  by 

®B  “  II 

where  X  is  the  optical  wavelength  and  A  is  the  acoustic  wavelength. 

This  type  of  interaction  resuits  in  only  a  single  first  order  diffracted 
beam  and  the  portion  of  the  origianl  beam  which  is  not  diffracted.  In 
the  figure  shown,  the  effect  of  refraction  according  to  Snell's  law  is 
omitted  for  simplicity. 
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The  interaction  between  the  light  and  acoustic  waves  may  fall 
into  one  of  several  different  regimes  depending  upon  the  optical  and 
acoustic  properties  of  the  material  used,  the  optical  and  acoustic 
frequencies  used,  and  parameters  of  the  cell  and  input  transducer*  in 
the  Raman-Nath  regime,  also  known  as  the  Debye-Sears  or  Lucas-Bequard 
region,  the  light  enters  parallel  to  the  acoustic  wave  fronts  and  many 
positive  and  negative  orders  of  diffraction  are  produced.  The  Raman-Nath 
type  of  operations  are  limited  to  low  frequencies  and  small  bandwidth. 
(Ref.  16:131 )•  Increased  bandwidth  and  higher  frequency  operation  re¬ 
quirements  have  led  to  widespread  use  of  the  Bragg  regime.  This  type 
of  operation  which  was  shown  in  Fig.  8  is  characterized  by  the  require¬ 
ment  that  light  must  enter  at  the  Bragg  angle  and  only  a  single  diffrac¬ 
ted  beam  is  produced.  Recently,  a  new  group  of  materials  for  which  other 
regimes  of  operation  are  possible  have  been  investigated  for  use  in  OSP. 
Birefringent  materials  exhibit  different  indices  of  refraction  for 
differently  polarized  beams  of  light  and  the  properties  differ  according 
to  different  directions  of  propagation  relative  to  the  crystal  struc¬ 
ture.  Birefringent  crystals  present  greater  versatility  and  improved 
characteristics  for  signal  processing.  The  birefringent  cells  are 
capable  of  high  frequencies  and  large  bandwidths.  Either  a  single  first 
order  or  positive  and  negative  first  order  outputs  may  be  produced 
depending  on  the  construction  of  the  cell. 

The  essential  properties  of  the  interaction  taking  place  in  all 
acousto-optic  cells  may  be  explained  by  considering  the  interaction  as 
a  collision  between  photons  and  phonons  (Ref.  21:49).  Conservation  of 
momentum  is  satisfied.  The  moments  of  the  particles  are  given  by 
11  k  and  fi  K  for  the  photons  and  phonons  respectively,  where fi  ■  h/2iT 


and  h  is  Planck's  constant.  The  propagation  vectors  of  light  and  sound 
and  l<,  have  magnitudes  2TrA  and  2it/A  ,  ,  respectively  and  their  direc¬ 
tions  are  the  directions  of  propagation.  From  the  wave  vector  diagrams 
in  Fig.  9  for  Bragg  diffraction  the  Bragg  angle  is  seen  to  be  given  by 


I  ~  I  K  1  X 
sin  ttB  -  y  y  -  y  y 


as  given  earlier  in  (12).  Both  down  shifted  and  up  shifted  interactions 
are  shown.  These  shifts  can  be  found  from  the  conservation  of  the  photon 
and  phonon  energies  liai,  and  (Ref.  15:49).  The  angle  at  which  de¬ 
flection  takes  place  is  related  to  the  frequency  content  of  the  acoustic 
wave  by  the  magnitude  of  j<.  The  cell,  therefore,  naturally  analyzes 
the  spectrum  of  the  signal  input  to  the  cell.  Note  that  for  the  inter¬ 
action  to  take  place,  K,  and  thus  the  frequency  input,  must  be  within 
a  limited  range.  The  proportion  of  the  input  beam  which  is  deflected  is 
related  to  the  power  contained  in  the  acoustic  wave. 
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The  relative  tight  intensity  of  the  first  order  diffracted  tight  beam 
for  Bragg  diffraction  has  been  given  as  (Ref.  17:54) 

-p  «  C  sinc^  I  C  +  (Ak  L/2)^  |  *  (14) 

o 

where  1^  and  1^  are  the  intensities  of  the  diffracted  beam  and  the  input 
beam,  P^  is  the  acoustic  power,  Ak  is  the  momentum  mismatch  of  the 
tight  and  acousto-optic  propagation  vectors,  is  the  width  of  the 
acoustic  beam  and  C  is  constant  which  depends  on  the  proper  ties  of 
the  ceil  materiai  and  its  dimensions.  Obviously,  the  response  of  the 
ceil  is  not  in  general  linear  either  for  amplitude  or  intensity  modula¬ 
tion..  To  achieve  linear  operation,  the  ceil  is  operated  only  over 
limited  ranges  and  efforts  to  reduce  the  effects  of  the  momentum 
match  Ak. 

The  design  parameters  associated  with  AO  ceils  are  efficiency, 
bandwidth,  resolution,  distortion  and  dynamic  range.  The  efficiency 
of  the  ceil  is  measured  by  the  percentage  of  the  input  beam  which  is 
deflected  per  watt  of  acoustic  input.  The  input  transducer  is  generally 
limited  to  less  than  one  watt  so  a  ceil  which  has  poor  efficiency  can 
limit  the  optical  ouptut  for  a  given  input  light  source.  AO  ceil 
efficiencies  range  from  1  to  300  t/W.(Ref  10:2).  The  bandwidth  limita¬ 
tions  of  the  ceil  are  due  to  a  number  of  sources  including  the  input 
transducer,  (Ref  I6:i83)»  the  width  of  the  optical  beam  and  acoustic 
beam  characters! tics  resulting  in  momentum  mismatch.  AO  ceils  have 
been  fabricated  with  bandwidths  of  40  to  1000  Mhz.(Ref  13:2).  The 
techniques  used  to  Increase  bandwidth  are  varied,  however,  in  general 
bandwidth  increases  can  only  be  obtained  by  sacrificing  other  performance 
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characteristics  such  as  aperture  and  efficiency. 

The  resolution  of  the  AO  cell  Is  Important  to  the  performance  of 
the  TIOC.  This  performance  Is  affected  by  the  Information  capacity  of 
the  cell  which  can  be  approximated  by  the  length  of  the  time  window 
which  is  the  aperture  multiplied  by  the  bandWi/dth  of  the  cell.  The 
resolution  is  further  limited  by  diffraction.  Elements  which  are  re¬ 
solvable  at  the  AO  ceil  any  not  be  resolvable  at  the  detector  due  to 
diffraction.  The  diffracted  image  has  reduced  resolution  due  to  spread¬ 
ing  of  the  beam  and  also  the  si delobes  produced.  The  effects  of  the  side- 
lobes  can  be  reduced  by  apodizatlon  of  the  optical  beam  (Ref.  16:136). 

The  resolution  of  the  cell  is  also  improved  by  increasing  the  aperture 
and  the  bandwidth  thereby  increasing  the  time  bandwidth  product,  however, 
as  mentioned,  increased  bandwidth  is  usually  accompanied  by  shortened 
apertures..  For  the  purpose  of  this  thesis  it  will  be  assumed»as  is 
the  common  practice  (Ref.  1:33)fthat  the  number  of  resolvable  points 
in  the  image  plane  is  equal  to  the  timeband  width  product.  AO  cells 
have  been  produced  with  timeband  width  products  of  as  high  as  1100  for 
a  bandwidth  of  1.05  GHz  or  1600  for  kO  MHz  bandwidth. 

A  number  of  nonlinear  characteristics  are  present  in  the  AO 
cell  which  can  contribute  to  distortion  both  in  the  amplitude  or  in¬ 
tensity  of  the  modulated  output  and  in  its  phase.  The  sources  of 
the  limitations  and  the  techniques  used  to  counteract  them  are  too 
numerous  and  complex  to  be  included  in  this  thesis. 

The  dynamic  range  of  the  AO  cell  is  limited  by  the  range  of 
the  acoustic  input  transducer,  and  by  the  regions  of  operation  over 
which  the  overall  response  is  linear.  This  parameter  is  important 
because  of  the  restriction  it  imposes  on  the  signals  which  can  be 


processed  using  these  devices.  Cells  have  been  produced  for  which  dy¬ 
namic  ranges  have  been  claimed  to  be  as  high  as  80  dB;  however,  60-65 
dB  is  a  more  realistic  limit  for  today's  technology  (Ref.  10:2). 

Sources  and  Detectors 

Sources  and  detectors  for  TIOC  will  be  discussed  here  together 
since  the  random  fluctuations  in  the  light  sources  will  most  easily  be 
described  in  terms  of  the  noise  which  results  upon  detection.  The 
sources  which  will  be  discussed  are  the  laser,  laser  diodes  (LD),  and 
light  emitting  diodes  (LEO).  The  detectors  required  for  the  TIOC  are 
produced  as  arrays  of  photo  diodes  whose  outputs  are  integrated  and 
then  read  out  at  the  end  of  the  detection  period. 

Sources.  The  nature  of  the  "noise”  present  In  light  from  different 
sources  can  be  described  by  the  variance  of  the  photo  electron  count 
during  the  observation  time  T^.  This  noise  is  commonly  referred  to  as 
shot  noise.  For  an  ideal  monochromatic  (zero  line  width)  source  the 
counting  statistics  are  independent  over  disjoint  observation  periods 
resulting  in  a  counting  process  which  is  Poisson.  The  variance  of  the 
count  Q  is 

aj  -  E  (Q)  (15 

where  E{Q}  is  the  expected  value  of  the  count  over  the  observation 
period.  If  a  laser  is  designed  carefully  to  produce  only  a  single 
temporal  mode  the  line  width  may  be  made  very  narrow,  therefore. 


is  a  good  estimate  of  noise  variance  exhibited.  Most  iasers,  however, 
are  designed  to  produce  several  temporal  modes  to  increase  output  power. 
These  modes  are  at  different  wavelengths  and  their  random  phase  varia¬ 
tions  introduce  beat  noise  also  referred  to  as  excess  photon  noise 
(Ref.  18:696,19:732,  and  20:16).  This  produces  a  Bose-Einstein 
counting  distribution  similar  to  that  for  thermal  photon  generations 
where  the  variance  is 

aj  -  E{Q)  +  E^CQ)  .  (16) 

This  expression  is  valid  where  the  coherence  time  remains  much  greater 
than  the  observation  time  and  the  coherence  area  is  much  greater  than 
the  detector  area.  The  excess  noise  term  can  be  very  large  for  large 
counts.  This  situation  should  obviously  be  avoided  for  communications 
or  optical  signal  processing  applications. 

it  is  noted  by  Pratt  (Ref.  20:l6)  that  phase  locking  the  modes 
of  the  laser  eliminates  the  random  nature  of  the  mode  beating  returning 
the  iaser  to  Poisson  statistics. 

The  laser  diode  (LD)  is  typically  less  coherent  than  other 
forms  of  the  laser  resulting  from  the  large  number  of  laser  lines  or 
temporai  modes  which  are  exited.  The  bandwidth  of  the  laser  diode 
output  may  be  on  the  order  of  500  Ghz  containing  perhaps  twenty  laser 
lines  (Ref.  21:96).  The  light  emitting  diode  (LED)  emits  over  a  band¬ 
width  on  the  order  of  20-30  THz.  The  noise  variance  for  these  sources 
is  given  by 


(17) 


2  \ 

E{Q}  +  e2{Q}  ® 


V>Tc» 


A»A 

W 


where  A  is  the  coherence  area*  A  is  the  detector  area.  T  is  the 
coherence  time*  and  is  the  observation  period.  The  coherence  time 
for  a  multimode  LD  is  about  29  pico  seconds  or  less  (Ref  21:98)  and  the 
coherence  time  for  the  LEO  is  on  the  order  of  pico  seconds  which  is  about 
the  same  as  from  a  thermal  source.  For  reasonable  observation  times* 
the  ratio  T^/T^  is  very  small  and  the  ratio  A/A^  which  represents  the 
number  of  spatial  modes  received  is  very  large.  For  the  LED  the  excess 
photon  noise  is  certainly  negligible  and  depending  on  its  coherence  this 
term  may  also  be  neglected  for  the  LD  source. 

An  important  point  to  note  here  is  that  in  general  incoherent 
systems  using  LED  sources  use  intensity  modulation  of  the  light  output. 
Light  detectors  respond  to  incident  power  so  this  method  produces  linear 
output  response  to  the  input  signal.  Coherent  systems  on  the  other  hand 
may  use  amplitude  modulation  of  the  light  field  thus  requiring 
heterodyne  detection  at  the  output.  Even  In  coherent  systems  using 
direct  detection  difficulties  may  be  encountered  due  to 
phase  errors  introduced  by  the  optics  or  the  modulation  process  resulting 
in  interference  problems  and  degradation  of  system  performance. 


Detectors.  (Ref.  22).  The  detectors  required  for  the  TIOC  must  be  able 
to  resolve  at  least  as  many  points  as  the  number  of  resolvable  points 
in  the  AO  celt  and  they  must  provide  an  integrated  output.  These 
requirements  are  generally  approached  using  CCD  or  CID  photo  detector 
arrays  (charge  coupled  devices  or  charge  inspection  devices).  This 
section  will  present  a  brief  discussion  of  each  of  the  noise  sources 
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present  in  these  devices  and  discuss  some  other  specifications  pertinent 
to  the  TIOC  performance. 

The  temporal  noise  sources  in  the  detector  array  include  thermal 
noise  (Johnson  noise)  associated  with  charging  capacitors  through  a 
res i stance »  shot  noise  due  to  leakage  current ,  the  shot  noise  of  bulk 

I 

traps  in  MOS-CCOs  and  the  photon  noise  which  was  discussed  in  association 
with  light  sources  (Ref.  22:108).  The  detector  arrays  perform  integration 
by  first  charging  a  capacitor  associated  with  each  array  element  to  a 
known  value.  The  detected  photo-electrons  discharge  the  capacitor  and 
the  integrated  signal  is  determined  by  the  difference  between  the  ori¬ 
ginal  charge  and  the  charge  remaining  after  detection.  Both  charging  and 
discharging  (read  out)  of  the  capacitor  introduce  noise.  The  variance 
of  this  noise  is  directly  related  to  the  temperature  of  operation  so  cool¬ 
ing  is  sometimes  used  to  improve  performance.  The  leakage  current  or  dark 
current  introduces  shot  noise  with  Poisson  statistics.  This  adds  to  the 
variance  of  the  output  and  is  also  detrimental  because  it  discharges  the 
preset  charge  In  the  readout  capacitance,  thus,  limiting  the  dynamic  range 
of  the  device.  Another  noise  source  is  bulk  trapping  noise.  The  contri¬ 
bution  from  this  source  is  very  small  except  for  cooled  arrays  with  long 
integration  times  (Ref  22:109)  and  its  details  will  not  be  discussed  here 
except  to  note  that  it  presents  a  nonPoisson  shot  noise. 

An  important  source  of  noise  for  all  detector  arrays  is  fixed 
pattern  or  spatial  noise  due  to  nonuniform  response  across  the  array. 

This  noise  is  due  to  a  large  number  of  sources  within  the  detector. 

Post  processing  can  be  used  to  compensate  for  this  noise  but  this 
approach  is  not  practical  for  the  number  of  elements  needed  for  TIOC. 


The  effects  of  all  of  these  noise  sources  is  to  reduce  the  dy¬ 


namic  range  of  the  device.  Dynamic  range  is  defined  here  to  be  the  ratio 

of  the  maximum  photo  electron  count  detectable  which  is  limited  by  the 

initial  charge  placed  on  the  capacitor  and  the  minimum  detectable  photon 

count  which  is  limited  by  the  initial  charge  placed  on  the  capacitor  and 

the  minimum  detectable  photon  count  which  is  limited  by  the  noise  count 
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due  to  all  of  the  noise  sources  discussed.  Dynamics  ranges  of  10'^  to  10 
(cooled)  are  possible  with  current  technology.  It  was  noted  in  Chapter  il 
that  this  dynamic  range  limits  the  gain  of  the  TIOC.  This  limit  can  be 
overcome,  however,  by  multiple  readout  during  the  integration  as  well  as 
other  methods  (Ref  22:108). 

Summary 

This  chapter  has  presented  some  of  the  operating  characteristics 
of  AO  modulators,  light  sources  and  detector  arrays  as  they  are  used  in 
time  integrating  optical  correlators.  The  limitations  of  these  devices 
have  been  presented  and  for  many  of  the  important  paraments  current 
device  capabilities  have  been  cited.  The  numerous  variations  on  the 
architecture  of  this  processor  which  are  possible  have  not  been  discussed. 
However,  it  is  important  to  note  that  time  integrating  correlators  may 
be  constructed  to  produce  either  real  or  complex  correlations.  For  com¬ 
plex  correlations,  the  processor  must  either  perform  four  real  corre¬ 
lations  using  quadrature  components  of  the  complex  signals  or  the  pro¬ 
cessor  must  be  capable  of  detecting  both  phase  and  amplitude  of  the 
correlation  product  by  using  some  form  of  heterodyne  or  interferometric 
detection.  The  interferometric  detection  can  be  performed  by  using  co- 


herent  light  and  a  local  oscillator  or  it  might  also  be  accomplished  in 
an  incoherent  system  as  described  by  Keliman  (Ref  23).  Another  aspect 
of  this  processor  which  has  not  been  discussed  is  the  bias  levels  re¬ 
required  both  at  the  source  and  in  the  AO  celt  to  represent  both  posi¬ 
tive  and  negative  quantities.  These  bias  levels  are  different  for  each 
architecture  and  thus  cannot  be  specified  in  general,  however,  their 
effect  on  system  performance  is  to  reduce  the  dynamic  range  of  the  corre¬ 
lator  output  thus  limiting  the  overall  processor  gain. 

This  thesis,  having  noted  the  limitations  of  the  components  of 
the  time  integrating  opticai  correlator  will  assume  that  the  correlator 
response  is  linear  while  keeping  in  mind  the  resolution  and  dynamic 
range  limitations  for  which  this  is  true. 
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V.  A  Time  Integrating  Optical  Correlator 
Adaptive  Array  Processor 

Introduction 

The  previous  chapters  have  presented  the  background  of  adaptive 
arrays  and  optical  signal  processing.  This  chapter  will  present  an  AA 
processor  which  uses  the  time  integrating  optical  correlator  to  provide 
estimates  of  the  directions  of  arrival  and  power  of  up  to  N  noise 
sources.  The  direction  estimates  are  then  used  to  generate  complex 
weights  such  that  the  resulting  antenna  pattern  has  nulls  of  some  pre¬ 
set  depth  in  the  desired  directions.  After  presenting  the  basic  form 
of  the  processor  this  chapter  will  compare  a  maximum  likelihood  estima¬ 
tion  of  source  directions  with  this  processor  and  present  a  bound 
on  the  performance  of  the  system.  Also  included  in  this  chapter  are 
discussions  of  spatial  anguiar  ambiguity  and  resolution,  an  approach 
to  efficient  post  processing,  and  some  variations  on  the  form  of  the 
TIOC  for  AA  processing. 

The  Processor 

The  basic  ideas  of  the  TIOC  for  AA  processing  were  first  in¬ 
vestigated  by  Casasent  (Ref  1).  The  processor  forms  the  correlation 
of  the  signal  at  the  reference  element  (the  far  right  element)  with 
the  signal  at  each  of  the  other  elements  to  form  N  correlations.  For 
the  case  of  a  single  interference  source  the  signal  at  each  of  the 


elements  is  simply  a  time  delayed  version  of  the  signal  at  the  reference 
element  with  the  time  delay  given  by  (18) 


T  -  7  sin  0  (18) 

c 

From  Table  1»  the  signal  at  the  reference  element  from  interference 
source  k  is  and  the  resulting  signal  from  K  sources  is 

K 

ro(t)  -  Z  f^(t)  (19) 

The  delay  for  interference  source  k  at  an  angle  6|^  from  bores ight  at 
element  n  of  the  linear  array  with  evenly  spaced  elements  if  nT|^.  The 
signal  at  element  n,  for  n  »  1,2,.,.,N, 

K 

r_(t)  -  E  f^(t-n  T.  )  (20) 

n  k-1  ^ 


The  adaptive  array  processor  forms  the  N  correlations 


Each  of  the  correlations  formed  contains  the  following  terms 
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All  of  the  cross  terms  (k  ]<  I)  are  zero  for  sources  which  are 
either  independent  or  not  coherent.  The  remaining  terms  are  each  the 
correlation  of  a  source  signal  with  a  time  delayed  version  of  that 
signal.  These  correlations  are  combined  to  produce  N(d)  as  follows 

N  . 

N(0)  -  Z  R  sine)  (23) 

n=*l  "  ^ 

The  locations  of  the  peaks  in  N(e)  indicate  the  directions  of  the  inter- 
ference  sources  and  the  amplitude  of  each  is  proportional  to  the  power 
of  the  corresponding  noise  source.  The  TIOC  AA  processor  implements  this 
approach  by  performing  the  required  correlations  with  acousto-optic 
time  irftegrating  optical  correlators.  The  structure  of  the  processor  is 
shown  schematically  in -Fig.  10. 

Also  showr>  in  the  figure  are  the  output  correlations  which  would 
be  formed  for  a  single  Jammer  with  a  rectangular  waveform.  Note  that 
both  positive  and  negative  correlation  delays  can  be  produced  by  proper 
timing  of  the  reference  input  r^(t).  The  correlation  peaks  line  up 
indicating  the  direction  of  the  source.  The  discussion  so  far  has  con¬ 
sidered  only  directional  noise  sources  and  no  noise  in  array  elements  or 
in  the  processor.  In  the  next  section,  noise  will  be  added  to  the  con¬ 
sideration  and  signal  characteristics  will  be  discussed. 

Maximum  Likelihood  Comparison 

This  section  will  show  that  for  certain  conditions,  the  pro¬ 
posed  processor  provides  maximum  likelihood  estimates  of  the  source 
directions.  The  maximum  likelihood  estimate  of  a  vector  of  source 
directions  will  be  derived  and  this  will  be  shown  to  be  equivalent  to 
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Fig  10.  TIOC  AA  Processor  Structure  and  Output 


the  processor  output.  The  directional  Interference  sources  for  this 
derivation  will  be  assumed  to  be  deterministic.  In  addition,  the 
sources  will  be  teated  as  real  signals.  It  is  important  to  note  here 
that  this  processor  could  also  be  implemen-ted  to  operate  on  complex 
signals  by  performing  complex  correlations.  The  signal  from  the  k-th 
source  has  been  given  as  f|^(t).  The  signals  from  which  the  estimates 
are  to  be  obtained  will  be  taken  as  the  signals  at  the  elements  of  the 
array  S^(t),  corrupted  by  zero  mean  gaussian  noise  sources  n^^Ct)  with 
variance  where  n^(t)  is  the  noise  associated  with  the  n-th  element. 

The  noise  sources  are  assumed  to  be  independent  which  implies  that 
they  do  not  derive  from  directional  noise  sources.  As  can  be  seen  from 
Fig.  11,  which  shows  two  sources  and  a  three  element  array,  the  signal 
r^(t)  from  the  n-th  element  has  the  form 

■■.<')  ■  S.,  •’kl*  -  +  "n<'> 

If  '■-(t)  Is  represented  by  T  time  samples  where  r  ,  •  r  (t.)  and  the  out- 
puts  from  all  of  the  elements  are  combined  to  form  the  vector  £  which 
has  the  form 


Fig.  11.  Array  Signal  Structure 
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for  an  array  of  N  4-  1  elements.  The  desired  directions  can  be  derived 
directly  from  the  associated  time  delays  so  we  will  form  an  estimate  for 
the  vector  ^  whose  elements  are  the  K  element  to  element  time  delays 
T|^.  The  condition  probability  density  of  £  given  ^  is 


P(l/l)  -  TT 

n-0 
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The  maximum  likelihood  estimate  of  t  Is  the  vector  T  which  maximizes  the 
likelihood  function  A(t)  ■  In  (p  (£/t)).  This  is  equivalent  to  minimi¬ 
zing  the  magnitude  of  the  exponent. 

After  dropping  the  denominator  which  has  no  dependence  on  Xj,  the 
quantity  to  be  minimized  Is 
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After  performing  the  square  the  resulting  terms  are 
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The  first  term  does  not  depend  on  ^  so  it  will  not  affect  the  minimization. 
Terms  (i^a  and  29b)  can  be  rewritten  with  the  summation  over  time  replaced 
by  an  integral. 
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(30.) 
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At  this  points  it  is  alisumed  that  the  integral  of  the  cross  pro<ducts  in 
(30  )  is  not  dependent  on  T.  This  condition  is  generally  satisfied  if 
the  source*)  are  not  coherent  with  each  other.  Given  this  assumption  and 
noting  that  the  integral  of  the  squared  tterms  does  not  depend  on  ;r»  it 
is  observed  that  maximizing  the  likelihood  function  can  be  accomplished 
by  maximizing  the  magnitude  of  (30a).  Noting  that  the  n  equals  zero  term 
does  not  depend  on  j.  The  estimate  is  defined  as  ^  such  that 


N  /•  *  K 

n-i  "  k-1  ^ 


is  maximized.  This  can  be  compared  to  the  correlations  formed  by  the 
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proposed  processor  as  was  given  in  (22) 


if  r^  is  replaced  by  (2U) 
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Replacing  Tq  by  n  Tq  to  scale  each  correlation  as  indicated  in  (23)  and 
summing  over  n  the  TIOC  produces 


N(0) 
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where  the  dependence  of  Tq  on 6  is  given  by  (19)* 

If  the  noise  n^(t)  is  uncorrelated  with  the  directional  noise 
sources  (or  assuming  high  signa1-to*noise)  it  is  clear  that  determining 
the  values  of  Tq  which  maximize  N(9)  is  equivalent  to  maximizing  the 
likelihood  function.  It  is  to  be  noted  that  maximum  likelihood  requires 
that  (33)  be  maximized  independently  for  each  source  k  and  that  N(8)  sets 
all  T|^  values  equal  to  and  find  k  values  which  maximize  N.  The  opera¬ 
tions  are,  in  fact,  equivalent  for  wideband  or  incoherent  sources.  If 
the  sources  are  replaced  by  random  sources  derivation  of  the  maximum 
likelihood  estimate  is  mathematically  difficult;  however,  if  the  inte¬ 
gration  is  long  enough  the  correlation  operation  produces  a  good 
estimate  of  the  covarience  and  it  is  apparent  that  for  wideband  sources 
the  peaks  produced  in  N(e)  provide  estimates  of  the  directions  of  the 


sources 


Cramer-Rao  Bound 


To  quantify  the  performance  of  the  adaptive  array  processor,  it 
is  helpful  to  investigate  the  Cramer-Rao  bound  on  the  variance  of  the 
estimate  produced  by  maximum  likelihood  estimation.  The  bound  for  the 
real  signal  case  which  has  been  assumed  to  this  point  is  derived  here. 
The  derivation  for  complex  signals  would  be  similar  to  that  given  by 
Whalen  (Ref  24:337-3^2)  for  the  variance  of  a  time  of  arrival  estimator 
for  radar  signals.  The  Cramer-Rao  bound  is  given  by 

■  a-ln(p(r/T))  1.* 

{var  T/t}  >  E{^ - - - }  ^  (35) 

””  “L  31^ 

It  was  shown  previously  that  the  log  likelihood  function  was  proportional 
to 


remembering  that  the  derivative  Is  evaluated  at  the  true  value  of  x.  If 
f(t)  has  a  Fourier  transform  given  by  F(ju))  then 
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and  by  Parseval's  theorem 
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assuming  that  nearly  all  of  the  energy  In  f|^(t-nT|^)  Is  contained  In  the 
Interval  from  0  to  T,  The  derivative  Is  then  easily  evaluated  given  that 
the  Fourier  transform  of  3^(x(t))/3T^  Is  -  X(ja)).  The  resulting  bound 

on  the  variance  Is 
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The  energy  In  the  kth  signal  Is 

CO 

h-h  f  l»k(»l^d«>  ("') 

«CO 

The  noise  variance  can  be  replaced  by  N^/2  and  substituting  (4l)  Into  (40) 
yields 

[N  K  2E.  -  1-1 

var  {T/t>>|e  |  (42) 


where  0|^  is  a  measure  of  bandwidth  of  the  k-th  source  defined  by 

M(ju)  “  dw 
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as  given  by  Whalen  (Ref  24:339)*  This  bandwidth  measure  is  the  second 
moment  of  the  signal  spectrum  about  ui*  0.  Whalen  shows  that  a  similar 
problem  with  a  narrowband  signal  with  random  phase  yields  similar  results 
with  B  redefined  as  the  second  moment  of  the  signal  spectrum  about  the 
mean  of  the  spectrum  (Ref  24:342). 

The  important  conclusions  to  be  drawn  from  the  bound  as  expressed 
in  (42)  are  the  impr^evements  in  the  estimate  to  be  gained  from  higher 
signal-to-noise  and  larger  numbers  of  correlations  N.  Perhaps  even  more 
important,  is  the  observation  that  the  performance  of  the  estimator  and 
thus  the  processor,  under  consideration  is  dependent  upon  large  bandwidth 
signals. 


Amb i qu I ty 

This  section  analyzes  the  ambiguity  present  in  the  output  of  the 
TIOC  AA  porcessor.  The  ambiguity  in  detection  of  the  angle  at  which 
an  interference  source  is  located  is  dependent  on  the  parameters  of  the 
AA  and  on  the  characteristics  of  the  interference  source.  The  symmetry 
of  the  linear  array  limits  the  range  of  angles  to  90  degrees  either  side 
of  boresight  (a  source  detected  at  37  degrees  may  be  at  37  or  143  degrees) 
To  understand  the  other  limitations,  it  is  easiest  to  first  look  at  the 
ambiguity  in  each  of  the  correlations  and  then  translate  this  back  to 
N(6).  For  a  completely  incoherent  source  (this  implies  infinite  band* 
width)  only  one  correlation  peak  is  produced  no  matter  what  range  over 


which  the  correlation  delay  is  varied.  Host  sources  of  radiation  at 
frequencies  less  than  several  gigahertz  are  highly  coherent  so  their 
autocorrelation  functions  are  periodic.  A  source  which  is  periodic 


with  period  and  has  a  coherence  time  greater  than  the  correlation 
time  will  have  an  autocorrelation  which  is  periodic  on  the  delay  Tq. 

If  the  distance  between  the  elements  for  which  the  correlation  is  pro¬ 
duced  is  greater  than  the  wavelength  X  >  cT^^  of  the  source  then  more 
than  one  correlation  peak  will  be  produced  in  the  range  over  Tq  which 
corresponds  to  the  otherwise  unambiguous  range  of  6  from  -90  to  90 
degrees.  This  is  illustrated  in  Fig.  12  for  an  array  of  five  elements 
spaced  X/2  apart.  The  source  for  this  example  is  a  pulse  train  with 
period  T^.  The  dashed  line  indicates  the  delays  associated  with  90 
degrees  and  the  correlations  are  shown  for  a  source  at  90  degrees. 

There  is  no  ambiguity  beyond  that  caused  by  symmetry  for  but  as 

can  be  seen  in  the  figure, numerous  peaks  are  present  between  -90  degrees 
and  90  degrees  for  the  other  correldtTons.  For  no  ambiguity  to  be  added 
by  the  periodicity  of  the  correlation  the  separation  between  the  reference 
element  and  the  element  with  which  the  correlation  is  performed  must  be 
less  than  one  half  of  a  wavelength.  The  maximum  unambiguous  range  for 
a  source  at  a  given  angle  is  not  the  same  for  sources  at  different  angles. 
The  maximum  unambiguous  range  in  Is  ^  and  (1)  can  be  used  to  find 
the  corresponding  range  of  angles.  This  range  which  corresponds  to 

“  -  •In"'  -  sin*'  (t-  Tj))  (M) 

The  effect  of  this  ambiguity  is  reduced  in  the  TIOC  for  AA  by  the 
addition  of  several  correlations  together.  The  locations  of  the 


ambiguous  peaks  do  not  line  up  for  angles  other  than  the  correct  peaks. 

Resolution 

The  resolution  of  this  system  is  limited  by  the  impiementation 
with  optical  processing  and  by  the  characteristics  of  the  signais  being 
detected.  The  limit  imposed  by  the  implementation  is  the  finite  resolu* 
tion  of  the  correlation  at  the  detectors  as  discussed  in  Chapter  IV. 

With  existing  devices  this  limitation  restricts  the  range  of  applications. 
Casasent  discusses  this  limitation  as  applied  to  adaptive  phased  array 
radar  (Ref  1:33)<  The  resolution  required  in  the  AO  cells  which  have  set 
time  windows  of  .5  to  kO  ysec  is  determined  by  the  velocity  of  propagation 
of  the  interference  signal.  For  an  AO  cell  with  a  time  window  of  1  ysec 
and  a  time  bandwidth  product  of  1000  the  smallest  resolvable  time  delay 
is  1  ysec  To  provide  100  resolvable  locations  from  boresight  to  90 
degrees  would  require  an  element  spacing  of  .3  meters.  Casasent  suggests 
elimination  of  this  limitation  on  resolution  by  using  an  adjunct  antenna 
element  17.3  meters  from  the  reference  element.  This  solution,  however, 
could  lead  to  ambiguity  as  discussed  in  the  last  section  if  the  inter¬ 
ference  sources  are  coherent  over  a  period  of  time  comparable  to  the 
arrival  delay  of  the  signal  from  the  adjunct  element  to  the  reference 
element.  Addition  of  several  more  adjunct  elements  would  decrease  the 
effects  of  these  ambiguous  correlation  peaks  as  indicated  in  the  last 
section.  It  should  be  noted  here  that  simply  increasing  the  element  to 
element  spacing  for  the  linear  array  is  undersirable  because  of  resulting 
null  placement  limitations. 

The  other  factor  contributing  to  resolution  limitations  is  the 
width  of  the  correlation  peak  which  is  related  to  the  bandwidth  of  the 


ircrference  source  as  shown  in  Chapter  III.  For  simplicity,  the  inter¬ 
ference  source  once  again  will  be  assumed  to  be  a  pulse  of  duration 
and  bandwidth  B  *  1/2  T^.  The  correlation  peak  width  is  taken  as  *  i/B. 
The  relation  between  the  correlation  delay  and  angle  is  given  once  again 
by 


T  «  -  ^  sln0  (kS) 

The  effect  of  this  nonlinear  transformation  is  that  the  resolution  at 
different  angles  from  boresight  varies  greatly.  The  relationship  for 
the  width  of  the  peak  occurring  inN(6)  is 

A0  -  sin’ll  (T3  +  ^)}  -  sln’^  (t^  +  ^  )}  (46) 

For  an  element  separation  d  of  10  meters  and  a  pulsewldth  of  1  nsec  from 
a  source  at  10  degrees  the  angle  A0  is  1.75  degrees.  For  the  same  source 
at  90  degrees  the  angle  A0  is  19.9  degrees.  The  width  of  a  peak  at  angle 
0  is  given  by  the  approximate  relationship 


Post  Processing 

The  interference  source  directions  provided  by  this  processor 
must  be  efficiently  utilized  to  generate  the  complex  weight  values  which 
produce  nulls  in  those  directions.  The  antenna  response  B(0)  resulting 
from  a  given  set  of  weights  was  given  in  (3)  and  is  repeated  here 
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(48) 


'  B(0)  -  2  W  exp(j  -22^  slnS) 

n-1  ^ 

where  f  is  the  frequency  for  which  B(0)  is  valid.  Casasent  notes  that 
this  is  a  discrete  Fourier  transform  for  the  case  where  6  Is  replaced  by 
the  source  angles  6|^  and  f  by  the  frequencies  of  the  sources  f|^  (Ref  1:62). 
This  approach  to  computing  the  weight  vector  W  is  not,  however,  applicable 
if  the  sources  are  at  different  frequencies. 

Another  approach  to  determining  the  weight  vector  is  to  set  up 
a  system  of  linear  equations.  Each  equation  determines  the  response  of 
the  array  for  a  given  direction  at  a  particular  frequency.  For  an  N<f1 
element  array  N-t-f  equations  must  be  set  up  to  completely  determine  the 
weights.  A  matrix  £  is  formed  whose  elements  are 

^  ^,0  •  “«>(■<  ^  ’"•'m  (‘'9> 

I  • 

The  vector  matrix  equation 


B-PW  (50) 

must  then  be  solved.  To  maintain  a  desired  response  at  bores ight  to 
the  desired  signal  at  frequency  f^  the  first  row  of  P  is  defined  as 

P,  „  -  1  for  n-1,  N+1  (51) 

and  B^  is  defined  as  the  desired  response  at  bores ight.  The  inter¬ 
ference  source  directions  9^,  and  frequencies  f^,  are  used  to  complete 
P,  and  to  provide  nulls  at  these  angles  and  frequencies.  The  remaining 
elements  of  are  set  equal  to  zero.  The  resulting  vector  matrix 
equation  for  W 
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can  be  solved  by  finding  only  the  first  column  of  the  IP  matrix.  This  is 
true  because  the  vector  £  has  only  one  nonzero  element  which  is  its  first 
entry.  The  special  way  in  which  Is  formed  produces  an  array  with  a 
number  of  unique  properties.  It  seems  likely  that  if  these  properties 
can  be  put  to  use,  an  efficient  solution  of  the  vector  matrix  equation 
can  be  found.  Casasent  suggests  a  technique  for  solving  this  system 
of  equations  using  projection  mapping  techniques  (Ref  1:66). 

Variations 

This  section  will  briefly  present  some  varlationson  the  basic 
form  of  the  processor  which  extend  Its  capabilities  or  enhance  its 
implementation.  The  first  variation  to  be  discussed  is  suggested  by 
Casasent  to  generate  the  source  frequency  information  which  would  be 
needed  to  place  nulls  at  specific  frequencies  as  discussed  in  the  pre¬ 
vious  section.  The  processor  proposed  uses  time  and  space  integrating 
concepts  within  the  same  processor  to  produce  the  correlation  output 
already  discussed  and  produces  an  estimate  of  the  noise  power  contained 
different  frequency  bands  by  correlating  a  bandpass  filtered  signal 
with  the  original  signal.  The  filtering  operation  is  performed  using 
the  Fourier  transform  property  of  light  propagation  as  discussed 
in  Chapter  1 1 . 

One  other  variation  on  the  basic  form  of  the  processor  is 
shown  in  Fig.  13.  This  processor  uses  multiple  sources  with  a  single 
acousto-optic  cell  to  generate  the  desired  correlations.  The  sources 
are  at  different  wavelengths  Xj,  \2»  and  X^  and  are  modulated  by 
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the  signals  r|(t),  r2(t)  and  r^(t).  The  signal  from  the  reference  element 


r^(t),  is  used  to  modulate  each  of  these  and  a  prism  is  used  to  separate 


the  resulting  outputs.  The  products  r^(t-T)  rj(t),  r^(t-T)r2(t),  and 
r^(t-T)r^(t)  are  then  imaged  with  different  magnifications  onto  the  same 
detector  array.  The  output  then  has  the  form 


R(t) 


■A 


(t-Tj)r,(t)  +  r^(t-T2)r2(t) 


+  r^(t-Tj)rj(t)  dt 


If  the  lens  are  specified  correctly 


Tj  -  3  »  and 


./2  -  2  “^1 


(53) 


(54) 


This  processor  then  provides  a  sum  of  three  correlations  with  the  advan¬ 
tage  of  requiring  only  one  detector  array.  With  this  implementation,  the 
detection  noises  due  to  thermal  sources,  dark  noise,  quantization  error, 
and  readout  error  are  added  to  the  output  only  once  rather  than  3  times 
as  would  be  the  case  with  3  separate  correlators. 


VI.  Simulation  Results  and  Conclusions 


Introduction 

The  complex  nature  of  the  relationship  of  the  characteristics 
of  the  signal  environment  and  the  parameters  of  the  processor  to  the 
performance  of  the  system  prompted  the  development  of  a  computer  pro¬ 
gram  capable  of  simulating  the  processor  under  consideration.  The 
results  of  that  effort  are  presented  In  this  chapter.  The  simulation 
produced  results  which  confirm  the  conclusions  drawn  in  Chapter  V 
concerning  resolution  and  ambiguity.  Conclusions  concerning  a  com¬ 
parison  of  these  results  to  other  adaptive  array  processors  will 
then  be  given,  followed  by  some  suggestions  for  future  research 
in  this  area. 

Simulation  Results 

Simulation  of  the  TIOC  processor  and  its  environment  was  com¬ 
plicated  by  the  multitude  of  processor  and  signal  parameters.  The 
program  was  constructed  using  four  subroutines  called  by  a  main  program. 
A  listing  of  the  program  is  included  in  the  Appendix.  The  first  sub¬ 
routine  was  adapted  from  a  similar  subroutine  developed  by  Casasent 
(Ref  1:37)*  This  routine,  labeled  ARRAY,  generates  the  composite 
signal  present  at  each  element  of  an  arbitrary  size  linear  array  of 
equally  spaced  elements.  Any  number  of  Interference  sources  may  be 
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simulated.  Although  only  sinusoidal  and  periodic  rectangular  pulse 
signals  were  included  in  the  simulation  runs,  any  signal  source  for 
which  code  can  be  written  to  represent  a  function  of  time  could  be  sub¬ 
stituted  into  the  routine.  The  output  of  the  routine  is  a  matrix 
containing  a  vector  of  time  samples  representing  the  signal  at  each 
element  of  the  array. 

The  next  routine,  labeled  CORR,  performs  the  correlation  of 
the  signal  at  the  reference  element  with  the  signal  at  each  of  the  re¬ 
maining  elements  in  the  array.  The  output  is  a  matrix  containing  vectors 
of  correlation  values  at  different  correlation  delays  for  each  of  the 
correlations  formed. 

The  ANGL  routine  accepts  the  correlation  outputs  and  generates 
values  for  N(0)  at  every  one-half  degree  from  zero  to  ninety  degrees. 

These  values  are  then  analyzed  to  Identify  the  maximum  values  corres¬ 
ponding  to  the  angle  estimates  of  the  source  directions.  A  parameter  was 
Introduced  into  this  routine  called  RES.  To  help  prevent  identification 
of  several  points  associated  with  the  same  peak  only  those  points 
whose  angular  distances  from  previously  identified  points  are  greater 
than  the  self-imposed  resolution  limit,  RES,  are  retained  as  valid  points. 

The  final  routine,  WAITS,  accepts  the  angle  estimates  generated 
by  ANGL  and  determines  complex  weights  for  each  element  of  the  array 
by  solving  a  system  of  complex  linear  equations.  The  system  of 
equations  is  solved  to  generate  a  specified  response  at  bores ight  and  at 
each  of  the  estimated  source  angles.  An  IMSL  routine  LEQT1C  was  used  to 
provide  the  solutions  for  the  system  of  equations.  The  depth  of  the 
nulls  provided  by  the  computed  weights  Is  limited  only  by  the 
accuracy  of  LEQTIC.  WAITS  uses  the  computed  weights  to  determine  the 
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antenna  pattern  using  (3) 


Simulation  runs  were  performed  for  both  sinusoidal  and  periodic 
rectangular  pulse  interference  sources.  As  predicted  by  the  analysis 
in  Chapter  V,  the  resolution  for  the  sinusoidal  input  is  very  poor 
due  to  its  narrow  bandwidth.  The  processor  output  for  N(6)  with  one 
cosinusoidal  source  at  17  degrees  is  shown  in  Fig.  14  .  Poor  resolution 


Fig.  14.  Simulation  Output  for  One  Sinusoidal  input 


and  ambiguity  as  discussed  in  Chapter  V  degrade  the  processor  output 
for  the  sinusoidal  input  so  greatly  that  correct  angle  estimates  were 
produced  by  the  simulation  runs  for  only  the  single  source  case. 

This  case,  as  well  as  the  other  discussed  here,  were  performed  for  a 
five  element  linear  array  with  a  one-half  wavelength  element  spacing. 
Ninety  sample  points  were  produced  for  each  period  for  a  total  of  lOOi 
points  at  each  element.  The  integration  time  was  set  at  360  points 
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corresponding  to  four  periods  of  the  Input  waveform.  The  range  of  cor¬ 
relation  delays  was  set  at  400  or  Just  over  four  periods.  In  all  cases, 
the  period  of  the  interference  sources  were  equal  to  the  frequency  of 
operation  of  the  array. 

A  rectangular  period  pulse  train  signal  was  also  used  as  an  in  - 
terference  source.  This  model  for  the  source  allows  easy  variation  of 
the  bandwidth  by  varying  the  fraction  of  the  period  over  which  the 
pulse  is  present.  Shown  In  Fig.  15  sr®  the  results  obtained  for  a  source 
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Fig.  15.  Simulation  Output  for  Rectangular  pulse  Train  Inputs 

at  17  degrees  with  rectangular  pulses.  The  widest  peak  at  17  degrees 
corresponds  to  pulse  train  with  pulses  one-half  of  the  period  length 
which  will  be  referred  to  here  as  a  duty  of  0.5.  The  other  peaks  at 
17  degrees  correspond  to  duties  of  0.1  and  O.O5  for  the  intermediate 
and  narrowest  peaks,  respectively.  Note  that  these  sources  are  very 
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wideband  with  bandwidths  of  approximately  4,20,  and  40  times  the  operation 
frequency  of  the  array.  These  bandwidths  are  obviously  very  high  but  they 
are  illustrative  of  the  wide  bandwidths  needed  for  this  element  spacing 
to  achieve  good  resolution.  As  predicted  in  Chapter  V,  ambiguous  pulses 
are  also  present  at  52.0  and  72.5  degrees.  Due  to  the  fact  that  the 
results  from  four  correlations  were  combined,  the  false  peaks  are  lower 
than  the  correct  peaks.  Note  that  the  amplitude  of  each  of  the  sources 
used  for  this  example  was  adjusted  to  provide  the  same  energy  and  thus, 
about  the  same  peak  height  in  the  output. 

To  illustrate  the  reduced  resolution  at  angles  near  ninety  degrees, 
the  plot  shown  in  Fig.  16  was  produced.  This  is  the  simulation  output  for 
a  periodic  pulse  train  source  at  ninety  degrees  with  a  duty  of  0.1. 
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Wideband  inputs  allowed  the  simulation  program  to  be  run  for 
multiple  source  angles.  Fig.  17  shows  the  processor  output  for  four 


period  pulse  train  sources*  each  with  a  duty  of  0.2  at  angles  of  17*  45* 
53*  and  61  degrees.  The  peak  detection  routine  identified  peaks  at 
17*  45*  53*  and  62  degrees  and  the  WAITS  routine  generated  the  weights 
in  Fig.  17b.  it  is  to  be  noted  here  that  the  depth  of  the  nulls  in  the 
antenna  response  could  be  specified  at  any  value  limited  only  by  the 
hardware  used  to  implement  the  procedure. 

Simulation  of  a  great  number  of  other  cases  was  also  performed* 
however*  the  plots  which  have  been  presented  here  illustrate  all  of  the 
significant  results.  These  results  included  confirming  the  increase  in 
resolution  with  wideband  signals*  the  differences  in  resolution  at 
different  angles*  the  presence  of  ambiguous  peaks  for  periodic  sources* 
and  the  ability  of  the  processor  to  recognize  multiple  sources  and  to 
produce  the  proper  weights  to  null  these  sources* 

Processor  Performance  Evaluation 

Based  upon  the  adaptive  array  background  as  discussed  in  Chapter 
22*  the  TIOC  for  AA  processing  approach  can  be  evaluated  using  the 
results  of  the  investigation  of  time  integrating  optical  correlators  and 
the  analysis  of  the  characteristics  of  this  type  of  processor  as  pre¬ 
sented  in  Chapter  IV  and  V.  First  of  all*  it  is  important  to  note  that 
the  processing  approach  which  has  been  discussed  Is  a  higher  order  method 
in  contrast  to  mose  of  the  algorithm  which  uses  a  steepest  descent 
approach  (first  order).  This  allows  the  processor  to  adapt  in  a  single 
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processing  cycle  rather  than  approaching  the  optimum  In  small  steps. 
Another  point  in  favor  of  this  system  is  its  flexibility.  As  'scussed, 
the  processor  assumes  a  known  direction  of  arrival »  however,  the  system 
could  easily  be  modified  to  acquire  the  signal  based  on  correlation 
with  a  signal  having  the  desired  signal  characteristics.  The  processor 
could  obviously  also  be  operated  in  an  all  signal  suppression  mode. 

As  discussed  briefly  In  Chapter  V,  addition  of  spectrum  analysis  to  the 
system  could  provide  the  added  flexibility  of  producing  nulls  at  any 
specified  frequencies.  The  depth  of  the  nulls,  as  well  as  the  width 
or  frequency  range,  is  dependent  only  on  the  hardware  used  to  implement 
the  weights. 

The  problems  associated  with  this  processor  are,  however, 
numerous.  It  has  been  shown  in  Chapter  V  that  the  processor  is  the 
optimum  estimator  of  source  angle.  The  conditions  imposed  on  the  signal 
enviornment,  however,  limit  the  applicability  of  these  results.  These 
conditions  were  that  the  interference  sources  were  incoherent  over  the 
integration  time  thus  producing  n^xiperiodic  correlations  and  that  multi> 
pie  sources  were  uncorrelated  with  each  other.  As  shown  by  the  Cramer- 
Rao  bound  and  the  discussion  on  resolution,  performance  of  the  system 
is  also  highly  dependent  on  large  bandwidth  signals.  The  effects  of 
periodic  signals  which  are  coherent  over  the  integration  time  was 
discussed  In  the  section  on  ambiguity  and  the  ambiguous  peaks  were 
observed  in  the  simulations.  The  poor  resolution  created  by  narrow 
bandwidth  signals  was  shown  by  simulating  a  sinusoidal  input.  The 
effects  of  correlated  sources  have  not  been  analyzed  explicitly  here. 
However,  depending  upon  the  degree  to  which  the  sources  are  correlated. 
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it  is  obvious  that  additionai  peaks  would  result  in  the  correlation  out¬ 
put.  The  plot  in  Fig.  17  shows  many  extra  peaks,  some  of  which  are  due 
to  the  periodicity  of  the  correlations  and  some  ow  which  are  due  to  cross¬ 
correlations  between  the  sources. 

In  addition  to  these  limitations  imposed  by  the  array  environ¬ 
ment.  implementation  of  the  processor  with  time- integrating  optical 
correlators  sets  further  limits  on  resolution.  Large  array  element 
separations  are  required  to  overcome  this  resolution  limit.  The 
dynamic  range  of  the  processor  is  also  limited  by  the  numerous  noise 
sources  and  bias  levels  included  in  the  correlator  output. 

Taking  into  consideration  all  of  these  factors.,  it  appears 
that  for  certain  applications  the  TIOC  for  AA  processing  offers 
possibilities  for  flexibility  and  fast  accurate  performance.  The 
limitations  discussed  suggest,  however,  that  without  improvements  in 
the  resolution  of  the  time  integrating  optical  correlator  the  range 
of  applications  is  very  limited. 

Future  investigation 

There  are  many  possibilities  for  further  investigations  related 
to  the  TIOC  for  AA  processing.  Simulations  of  different  types  of 
sources  as  well  as  the  addition  of  nondirectional  noise  could  be 
added  and  various  element  spacings  could  be  simulated.  The  similarity 
of  this  processor  to  radar  lies  in  the  problem  of  detecting  time  of 
arrival  by  correlation.  Investigation  of  signal  processing  techniques 
used  to  enhance  radar  performance  should  be  made  to  determine  applicability 
to  this  processor.  Some  form  of  angle  gating  (corresponding  to  range 
gating  for  radar)  might  be  possible  to  Improve  performance  of  the  system. 


Summary 

This  thesis  has  investigated  the  application  of  optical  signai 
processing  to  adaptive  array  processing.  Adaptive  array  processing 
was  investigated  and  the  repeated  use  of  correlation  in  the  common 
algorithms  led  to  the  decision  to  concentrate  the  study  of  optical 
signal  processing  in  the  area  of  optical  correlators.  The  character¬ 
istics  of  both  spatial  integrating  and  time  integrating  correlators 
were  analyzed  and  sei action  of  the  time  integrating  form  was  made  based 
upon  the  signal  processing  needs  of  the  adaptive  array  and  the  capa¬ 
bilities  of  each  type  of  processor.  An  adaptive  array  processor  which 
was  first  investigated  by  Casasent  (Ref.  J  and  2)  was  presented  and 
this  processor  was  compared  to  an  optimum  processor  for  a  given  set  of 
conditions.  The  Cramer-Rao  bound  on  the  performance  of  this  adaptive 
array  processor  was  derived.  Ambiguity  in  the  output  and  the  resolu¬ 
tion  of  the  system  were  discussed  and  an  approach  to  post  processing 
was  presented.  Many  variations  on  this  processor  both  in  its  approach 
and  its  impiementat ion  are  possible  and  two  variations  were  discussed 
in  Chapter  V.  Simuiation  of  the  processor  and  its  environment  was 
accomplished  for  a  number  of  variations  and  the  resuits  for  severai 
cases  have  been  inciuded  in  this  chapter.  An  evaiuation  of  i imitations 
on  the  performance  of  the  time  integrating  optical  correlator  for 
adaptive  array  processing  for  different  applications  was  discussed  and 
it  was  concluded  that  in  many  situations  its  performance  would  not  be 
satisfactory.  The  positive  points  in  this  evaluation  suggest,  however, 
that  for  certain  enviconments  the  system  would  perform  well  and 


several  areas  for  future  study  of  this  processor  have  been  indicated. 
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APPENDIX:  SIMULATION  PROGRAM 


P RO 6R  AH  MAI N i INPUT t  OUTPUT tT APE9  > 

THE  MAIN  PPC6RAM  IS  A  UTILITY  TO  INITIALIZE  VALUES  AND 
TO  CALL  THE  SUBPROGRAMS 

DIMENSION  XR(10C0»5>«XSOURC20>tXFREQC20) 

DIMENSION  C0RPC400t5> 

REAL  MAN6L<20t2}tT(iq2) 

COMPLEX  QCStS) t6(3Gl>fWAITC5>«B(5> 

COMMON  /O/  NELEtPIfAOtSAMP 
COMMON  /A/  XR«NPTS«NSRC*BUMCD 
COMMON  /B/  TtCORRtNTAUtNCPTStNPEF 
COMMON  /C/  MANGLfQf GtUAITtBtFREQ 
COMMON  /E/  XSOURtXFREQ 
DATA  NELE«NPTStNSRC/5tlOOO»4/ 

DATA  XR/5000<0.01/*  FREQ/1.0/ 

DATA  XFRCa/20(1.0)/«XSOUR/17.0»45.0t53.Ct61.0fG7.9/ 

NREF  s  I 

NTAU=4()0 

NCPTS=360 

S<1 |:C10.0t0.0> 

B<2  }::B<3)>B<4)=PC5)=C0.0«0.0  ) 

DATA  AO/Q.5/ 

PRINT*t*NTAU  IS»tNTAU 
PRINT*»«NCPTS  IS**NCPTS 

PRlNT*f 'SOURCE  DIRECTIONS  ARE' t  CXSOURCI  >  tl  =  l  tNSRO 
PRI NT*f 'THERE  AFE  ' *NSRCf'30URCES' 

PRI Nr*f 'ANIENNA  HAS  'tNELE* 'ELEMENTS' 

PRI NT«f 'ELEMENT  SPACING  IS  • t AO ,' WAVELENGTHS' 

CALL  ARRAY 
CALL  CORREL 

PRlNT*f'ItJtCCRR(Iv J)'f <CI»J>CORR(It J)»I=1«400>» 

J=2«3I 
CALL  ANGL 

URI  TECtlSG  >(FLCAT(II/2.0  tl  (  I  )t  1=1 1 130  ) 
F0RMATClXt2HATflX»F7.3tlXfSH0EGREES=«E15.9> 

CALL  WAITS 
CALL  PLOTl 
CALL  PLOT 2 


SUaROUTINE  ARRAY 

SUBROUTINE  ARRAY  COMPUTES  NPTS  VALUES  FCR  THE  SIGNAL 
AT  EACH  ELEMENT  OF  THE  ARRAY  DUE  TO  THE  NSRC  JAMMING 
SOURCES.  THE  ELEMENT  TO  ELEMENT  DELAY  IS  FIRST  FOUND 
AND  THEN  THE  SIGNAL  AT  EACH  ELEMENT  AT  EACH  TIME  DUE 
TO  THAT  SOURCE  IS  COMPUTED.  THE  CONTRIBUTIONS  FROM 
ALL  OF  THE  SOURCES  ARE  SUMMED  AND  THE  RESULTS  ARE 
RETURNED  IN  THE  ARRAY  XR. 

DIMENSION  XSOURf20l«XFREQC20  >tXR(1000t5> 

COMMON  /□/  NELEtPItADvSAHP 
COMMON  /A/  XRfNPTSvNSRCtBUMOO 
COMMON  /£/  XSOURtXFREQ 
AMPL  =  0.1 

PI  =  3.1415926535898 
SAMP  =  90.0 
OUTY=0.1 
RUAVEL  =  l.G 

Pi»IMT*t*AMPLITUOE  IS  •*AMPL 

PRINT*f •THERE  ARE  • #SAMPf ’POINTS  PER  PERIOD* 

DO  100  IM  =  1«NSRC 
PHASE  =  2.0  •.PI  *  '»ANFf) 

PRINT*t*FOR  NSRC’f IM**PHASE  IS  ’tPHASE 
DO  90  I  =  1 »NPTS 
DO  BO  IL  s  IvNELE 
RNELE  =  FLOATCILI 

OELA><2.0*PI*AO*CRNELE-1.0>*SINC€f X30UR<IM>>»PII/Id0.0>>/RUAV 
OARG  =  <2.0*PI*CXFREQ«IM)>*r)/SAMP  ♦  PHASE  ♦DELA 
ISAMP=IFIXCSAMP) 

RDAR6=0ARG-FL0ATCIFIXf0ARG/(2*PI)))*2*PI 

XAR=0.0 

IF(ROARG.GT.DUTY*2*PI>GO  TO  75 
IF(RDARG.LT.O.OIPRINT*#*ERROR  IN  ARRAY* 

XAR=1.0 

CONTINUE 

XR(IflL)  =  XR<IfIL>  ♦  XAR 

CONTINUE 

CONTINUE 

CONTINUE 

URI  TEC*  f  120  MCXRf  It  J>*J=lt5>«I=lfl9} 

FORMAT(lXtlOE13.7) 


RETURN 

END 
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nan  n n n n n n n n n  nnn 


SUBROUTINE  CORREL 

SUBROUTINE  CORREL  PERFCMS  THE  CORRELATION  OF 
THE  SIGNAL  AT  ELEMENT  1  WITH  THE  SIGNAL  AT  EACH 
OF  THE  OTHER  ELEMENTS.  THE  INTEGRATION  IS 
PERFORMED  AS  THE  SUM  OVER  NCPTS  POINTS  OF  THE 
PRODUCT  OF  THE  SIGNALS.  THIS  IS  DCNE  WITH  A 
DELAY  ON  THE  SECOND  SIGNAL  CF  FROM  I  TO  NTAU 
POINTS. 

REAL  T(182) 

COMMON  /D/  NELEtPIt ADvSAMP 
COMMON  /A/  XR  tNPTStNSRC«BHMCD 
COMMON  /B/  TtCORRtNTAUtNCPTStNREF 
REAL  XR(1000t5)tCCRRC400«5> 

DO  60  I=ltNTAU 
DO  50  J^ltNELE 
CORRCItJI  =  0.0 
50  CONTINUE 
60  CONTINUE 


C 


1£0 

190 

200 


DO  200  lELE  =  ItNELE 

IF(  lELE  .EQ.  NREF  I  GO  TO  200 


1 

2 


DO  190  ITAU  =  IvNTAU 
DO  100  IT  =  1*NCPTS 

CORRCITAUtlELE)  =  CORRC I TAU * lELE )  ♦  .  XR (I T t NREF> •XR Cl !♦! TAUt lELE 

XRf IT*ITAU*IFIXCSAMP>»IELE>*XRCITtNREF)  ♦ 
XftCIT»ITAU*2*IFIX(SAMP) tIELE) *XRCIT*NREFI 


CONTINUE 

IFCITAU  .GT.  5>  GO  TO  190 

CONTINUE 

CONTINUE 

RETURN 

END 
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SUBROUTINE  ANGL 
C 

C  SUBROUTINE  ANGL  FINOS  THE  CCRRELATICN  AS  A  FUNCTI3N  OF  ANGLE 

C  AND  THEN  THE  SECOND  PART  OF  THE  ROUTINEt  THE  MANGL  ROUTINEt 

C  FINDS  THE  ANGLES  AT  WHICH  THIS  FUNCTION  TCTHETA)  IS  MAXIMUM. 

C 

COMMON  /O/  NELEtPI*AO*SAMP 
COMMON  /B/  TtCORR*NTAU.NCPTSfNREF 
COMMON  /C/  MANGL«QtGtUAIT>BtFREQ 
REAL  MANGL(20«2) 

COMPLEX  Q(5«5>tG(361)«WAITfS>*Bf5> 

DIMENSION  Tfl62)*CORRC400t51 #TC0PY(182> 

C  ELEMENT  1  IS  REF  ELEMENT 

C  ASSUMES  SOURCE  FREO  =1.0 

DO  200  I=ltl60 
TCI  »  =  0.0 

DELA  =  2.0*PI*AD*SINCFL0AT(I )»PI/360.01 
00  100  J=2«NELE 

TAU  =  (FLOATCJI-1.3 >*DELA*SAMP/(2.0*PI*FREQ) 

ITAU  =  IFIX<FL0AT<NCPTS)/SAMP>*IFIX(SAMP>  -  IFIXCTAU»0.5> 
IFCCITAU  .LE.  0)  .OR.  IITAU  .GT.  NTAU))  GO  TO  50 
TCI  I  =  TCri*CORR<irAU*J) 

60  TO  100 

50  URI TEC*«35}ITAUf Jf I 
55  FQRMATClXt5HERRCR«IS«2X*I3»2XtI3> 

100  CONTINUE 
200  CONTINUE 
C 

C  MANGL  ROUTINE 

C 

00  225  I=lt2 
DO  220  J=1«NELE 
MANGLCJtl 1=0.0 
220  CONTINUE 
225  CONTINUE 
C 

c 

DO  250  1=1*180 
TCOPYCIIsTCII 
250  CONTINUE 
RES  =  4 

DO  500  I  =  l*r.ELE 
300  CONTINUE 

MAN6LCItll=TCll 
MAN6LCI  *2>  =  C.5 
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00  400  JslvlBQ 

IF(r(J)  .67.  T<INTCMANGL{I»2)*2.0>)>  M AN6LCI »2 > -FLOAT CJ) /2.0 
400  CONTINUE 

T<INT(MAN6L(If2)*2.0)>  =  -10**20 
00  450  L=ltf  ELE 
IF(L.EQ.I>  GO  ^0  450 

IF( ABS(NANGL<It2>-NAN6L€Lt2l >  .GE.  RESIGO  TO  450 
GO  TO  300 
450  CONTINUE 

MANGLCIfl}=TCHANGLf I«2)} 

500  CONTINUE 

PRlNT*t«I  tHAN6L<I  t2l««(ItMA‘.6Lf  1 12)  tlsltNELE) 
PRINT*««I«MAN6L(Itt I • • < I t HANGLf I « 1 > • 1-1 t NELE) 

00  600  i^itiao 
T<I l=TCOPY<Il 
600  CONTINUE 
RET  URN 
END 


SUBROUTINE  WAITS 
C 

C  SUBROUTINE  WAITS  ACCEPTS  THE  MAXIMUM  ANGLE 
C  INFORMATION  FROM  THE  MANGL  ROUTINE  AND  USES 

C  THESE  ANGLES  TO  SET  UP  A  SYSTEM  OF  COMPLEX 

C  EQUATIONS  ^^PROOUCE  NULLS  AT  THE  FIRST 
C  (NELE-ll  ANGLES  LISTED  IN  MAN6L(I«J>«  A 

C  CONSTRAINT  IS  ACCEPTED  FROM  THE  MAIN  PROGRAM 

C  IN  THE  Bill  VECTOR.  THIS  CONSTRAINT  SETS 

C  THE  ANTENNA  RESPONSE  AT  BORESIGHT  TO  10.0  ♦  JO.O  • 

C  THE  ROUTINE  CALLS  AN  IMSL  ROUTINE  •LEOTIC* 

C  WHICH  SOLVES  THE  SYSTEM  OF  EQUATIONS  FOR  THE 

C  COMPLEX  WEIGHTS.  THESE  WAITS  ARE  THEN  USED 

C  TO  COMPUTE  THE  RESULTING  ANTENNA  RESPONSE  AT 

C  ANGLES  FROM  1  TO  360  DEGREES. 

C 

REAL  UA(5I«MANGL(20 t2) 

COMMON  /D/  NELEfPIt AOtSAMP 
COMMON  /C/  HANGLtQtGtWAITfBvFREQ 
COMPLEX  Q<5»5lfGC361)tUAITC5>»BC5)«X 
C 

DO  100  I=l«361 
GCII  =  (Q.OtO.Q) 

100  CONTINUE 
C 

DO  IDS  IslfNELE 
Qdtl)  =  fl.OtO.Ol 
C 

105  CONTINUE 


74 


DO  120  I  =  ltK‘ELE-l 

THETA  =  2.0*PI*FREQ<*SINfMANGL(It2>*PI/130-0) 

00  110  JsltNELE 

Q(  t^lfJ  )=CHPLXC  CGSC  FLOATC  J  ) *THET A» tSIN (  FL0AT<  J  )*THET 

CONTINUE 

CONTINUE 

PRiriT**»Q  MATRIX* 

PRINT*t(C*I  tJtQ<I««l>=  •vItJtQCI  tJ)tI-lt5>f  J=lf5) 

CALL  LEQTlC<QtNELEtNELE«Bfl «NELEtO»UAtIER) 

URITe(*tl3Q>CB(I)tI=ltNELEI 
FORMATIlXfllHUEIGHTS  AREt  fE15.9  tSKt E15 .9)  > 

COMPUTE  THE  ANTENNA  PATTERN 
GENERATED  BY  THE  COMPUTED  WEIGHTS 

DO  220  I=lf361 

THETA  =  2.0*PI*FREQ*SINCFLOAT(I-1)*PI/180.0> 

DO  210  JsltNELE 

X  ^MPLXCCQS  (FLOATC  JI*THETA)  f  SIN(FLOAT(  J>  *THETAI  > 

G(I)  =  6(11  ♦  B(J1*X 

CONTINUE 

CONTINUE 


DO  240  Islt361 
URITE(*«23Q)  X>lt6(I> 

F0RMAT(lX«2HATtI3«7HDE6REESt3X«£15.9*5XtE15.91 
CONTINUE 
RETURN 
END 
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20.  ABSTRACT  /Conllmia  on  roaoraa  alda  If  nacaaaofr  md  Iddntlty  by  block  numbarj 

Application  of  optical  signal  processing  techniques  to  adaptive  array 
processing  was  investigated.  The  common  adaptive  array  algorithms  were 
reviewed.  Adaptive  array  processing  generally  requires  some  form  of  correla* 
tion  and  the  characteristics  of  optical  correlators  were  compared  leading 
to  the  selection  of  the  time  integrating  optical  correlator  for  further 
investigation.  An  adaptive  array  processor  using  time  integrating  optical 
correlators,  which  was  first  proposed  by  David  Casasent,  was  analyzed  to 
determine  applicability  of  this  signal  processing  technique  to  a  general 
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adaptive  array  environment.  This  processor  was  found  to  be  similar  to  a 
comparable  maximum  likelihood  approach.  Limitations  on  the. performance  of 
this  processor  ia  terms  of  omblguity  and  resolution  were  found  and  these 
results  were  confirmed  by  computer  simulation  of  the  adaptive  array 
environment,  the  processor,  and  the  resulting  antenna  response. 
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